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ABSTRACT 


Data for planform (meandering) and depth profiles for three 
Canadian rivers (Beaver, Alberta; Kootenay, B.C.; South Nahanni, N.W.T.) 
have been collected in an attempt to better understand river behavior. 
Earlier workers attempted to explain such forms by secondary flows, 
transverse oscillations, and the earth's rotation, but did not predict 
river behavior well. More recently, stochastic models produced variance- 
minimizing sine-~generated curves which approximated meandering. Spectral 
analysis has been applied by Speight (1965) and others to analyze meanders, 
and by Nordin (1966) and other to analyze bedforms. Stochastic processes 
and variance spectrum analysis are reviewed. 

Collection techniques for the data acquired are outlined, as are 
the relevant physical geographic characteristics of the areas and rivers 
studied. Methods of digitizing, correcting, storing and analyzing the 
meander and depth series are presented. 

All series are at least weakly self-stationary. Digitizing and 
analytic techniques used are efficient and inexpensive. Simple statistics 
can reveal general trends in such data, but spectral analysis shows: 
these pronounced (down-path) wavelengths are between 6-20 river widths; 
most wavelengths occur as multiples of (#7 x width); and most spectra 
resemble those characteristics of a second-order autoregressive process. 

Recommendations for future research include: improvement in types 
of data collected; improvement in resolution of data actually used in 


analysis; better absolute location when collecting data; better data 
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South Nahanni River braided section 
(photo, chart, correction) 


"The Gate" in Third Canyon on South Nahanni River 
(photo, chart, correction) 


"Figure-eight" Rapids on the South Nahanni River 


Rapids on the South Nahanni River 


Transducer mount for the South Nahanni River 


Transducer in operation 


Standing Waves on the Kootenay River 


Research Council of Alberta jetboat 


Determination of 9 and R for Planform 
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Plan and Depth Spectra, Coherence 
upstream meanders, Beaver River 
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LIST AND EXPLANATION OF SYMBOLS IN THIS THESIS 


wave length 

lag specification 

total number of points in a series 

series mean 

cross-covariance function 

cross-correlation function 

cross-spectrum 

quadrature spectrum (imaginary component of the cross-spectrum) 


cospectrum (real component of the cross-spectrum) 


measure of the phase relations 


phase angle 

coherency spectrum 

sampling interval 

angular deviation along the flow path 

"relative curvature" - changes in flow direction assigned either 
positive or negative values to specify right vs. left turns 

in flow 

“absolute curvature" ~- absolute value of angular deviation - flow 
change considered only with respect to magnitude, ignoring 
specific direction 

radius of curvature 


river width, defined further in context 
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river depth, defined further in context 

data for South Nahanni River, upstream survey 
data for South Nahanni River, downstream survey 
data for Kootenay River, upstream survey 


data for Beaver River, upstream survey 
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INTRODUCTION 


A. The Problem 

Effective control of the behavior and form of natural systems 
can only be achieved through sufficient understanding of the various 
processes interacting to maintain the systems concerned. A system must 
be described, defined, and understood before it can be safely tamed or 
made more productive. Much remains to be learned about the river, and 
the potential contained in such knowledge is great. 

The tendency of all rivers to establish and maintain curved sections 
of channel in plan continues to fascinate fluviologists. The process of 
meandering is neither fully described nor fully explained. Even less data 
is available on the "meandering" pattern assumed in the vertical dimension 
by the river's bed. Laboratory work in fluvial geomorphology allows a 
variety of hypotheses to be tested against control situations. Nature 
does not offer the same advantages, and relatively little field data 
of great detail exists for natural rivers. Correspondence of lab and 
field studies must continually be tested toward the building of theory 
in fluvial geomorphology. 

Additional expansion of the potential utility of detailed field 
data, through the sophistication of methods in data collection and 
treatment, must be maintained as a major aim of work in fluvial geomorphology 
as theory is established. New techniques may be introduced from other 
fields or developed within the discipline, but the level of analysis of 


data must constantly be advanced. 
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Briefly, then, both the quality and quantity of data gathered 
towards solution of the problem of fluvial meandering must be improved. 
Analytic techniques must also be upgraded so that the level of under- 
standing is increased. Comparative utility of data and analyses must be 


maximized within the field. 


B. Purpose of Research 


Although it is possible to obtain continuous depth profiles of 
navigable rivers, previous studies in fluvial geomorphology have seldom 
exploited the potential of such data. Consequently, comparative studies 
using continuous depth profiles for rivers, covering a variety of channel 
types, bank and bed materials, geologic histories, river sizes, sediment 
loads, slopes, and discharges are undertaken for the first time at this 
scale in this thesis. Among the initial purposes established for this 


research were thus: 


1. to gather continuous depth profiles for rivers of different channel 
types, bank and bed materials, geologic history, size, sediment load, 
slope, and discharge, and 

2. to digitize, store, and so make available the data collected, for 
use by other researchers, government agencies, or any concerned 


individuals or groups. 


Comparative potential of this research was also considered 
important. Dr. M. Church of the University of British Columbia, is 
carrying out complimentary research, and he is primarily interested in 
scales of turbulence operating within fluvial systems. He has collected 


data on very small rivers(small mountain streams) and very large rivers 
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(the MacKenzie River). Data on rivers of intermediate size wiil thus 


complete the initial examinations of river plan and depth inter-relation 


oO 


with respect to scales of turbulence: 


3. to accumulate data on rivers of intermediate size for use in future 
research on plan and depth inter-relations for rivers of all sizes; 


GONUHEGLEN In conjunction with Ur. M. CGhtirch or Ura.G. 


Unless analysis of data can handle the detail gathered, its value 
remains. locked in potential utility. Superiority of analytic techniques 
with respect to quantity of data input and quality of information output 
is sought. More objective analysis must also be attained. Techniques 
employed in previous studies, where "typical" meanders were described, 
must be improved upon. Variance spectrum analysis provides sophistication 
not yet fully explored in conjunction with fluvial geomorphologic studies. 
Many problems are yet to be both discovered and subsequently solved in 
order for the technique to become established as a form of communication 
and adequate basis of comparative study in the field. Thus, other 


purposes underlying this research include: 


-. CoO atcempe Minimization of subjectivity, inherent in previous 
analyses of meander planform, and provide descriptions of the rivers 


StuldUueawOneehdise Dasass. 


5. to explore the utility of spectral and cross-spectral analytic 


techniques for data collected along the path of flow in rivers. 


6. to derive information from the data collected on both the planform 


geometry and bedform configuration of rivers, through the use of 
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spectral and cross-spectral as well as other "“less-sophisticated" 
techniques, and to explore the inter-relations of planform and 


depth, in order to provide descriptions of the rivers. studied. 


GC. scope Of the Thesis 

Availability of equipment, finances, and time constituted the 
major constraints encountered during research for this thesis. The use 
of depth recording machinery and boating equipment had to be arranged in 
co-operation with external organizations, and in this regard, research 
schedules became fixed rather than spontaneous. This limited examination 
of processes to those operating only at specified times, chosen rather 
arbitrarily, within the river systems studied. The study of processes 
operating at specific discharge stages (eg. bankfull) was thus beyond 
the scope of this study. 

Adequacy of first-hand knowledge of study areas also limited 
the scope of the study. Thalweg soundings were not recorded during the 
South Nahanni survey, as it was believed that this would create probleus 
with respect to safety of crew and equipment. It seems, in retrospect, 
that thalweg depths could actually have been obtained, but this knowledge 
was not available prior to the study. 

The scope of the field of research conducted for this study thus 
allowed only single samplings of the depth profiles of three accessible 
rivers. Lessons derived from the experience gained in the course of 
the research completed are listed in Chapter 5, under ''Recommendations 


for Future Research". 
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CHAPTER 1 


THEORETICAL CONTEXT OF THIS STUDY 


Early studies of the morphogenesis of fluvial landscapes 
immediately indicated the fluviologist's fascination with and underlying 
curiosity concerning the meandering form so consistently assumed by 
rivers. The physical behavior of the dynamic systems that mold the 
sinuosities so characteristic of fluvial sceneries is still not fully 
understood. Satisfactory statistical, geometrical, and descriptive 
definitions of the meandering form have not been attained in any wide 
sense within the discipline, but various approaches have each gained the 
acceptance and support of a number of fluvial geomorphologists during 
the development of the discipline. 

Classical explanations of the initiation, persistence, and form 
of fluvial meandering have varied widely in their approach to the problem. 
Examinations of form have been attempted in order that process might be 
understood, and process has been examined to explain form. As well, 
many physical properties have been examined simultaneously. Though many 
relationships have been established as hydraulic laws, and though many 
fluvial processes are now understood as complex components of river 
behavior, no single set of theories and no single model sufficiently 
explain the intricacies of meandering. 

A full appreciation of the theoretical problems and challenges 
involved in the attempt to explain the meandering tendency of rivers 
is best gained through an examination of the theoretical development 


to date. Many problems have been recognized in theory, and many have 
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1yition Of generating processes has 
been greatly increased and models have been substantially impeeved. 

Much of the subjectivity of the early analyses has been eliminated 
through sophistication of methodology and research technology. As 
theory has developed in other disciplines, fluvial geomorphology has 
grown, incorporating concepts of turbulence, and incorporating many 
ideas initiated in hydraulic engineering. Thus, although the meandering 
of rivers is not yet fully explained, considerable theoretical advances 
have been made. 

The principal theoreticians who have contributed insight into 
the problem of meandering are discussed below. In order to provide as 
objective an introduction as possible, contributions are presented in 
the terminology of the authors. Some criticism is then offered as the 


development is traced. 


A. Early Models 


Boussinesq (1872) offered an analogy between the flow in a bent 
pipe and that in a curved section of open channel. Thompson (1876) 
developed the physics of this concept and presented his own impression 
Gf flow, theorizing that, stretched out, a stream was°really screw- 
shaped so that flow was characterized by "transverse circulation". Though 
Thompson early atheveuau the phenomenon of transverse circulation, 
(later to become known as secondary flow), his theories did not achieve 
widespread support for some time. Nearly all early work in the fields 
of fluvial geomorphology and river control, extending into the early 
twentieth century, was based on the assumption that streamflow 
proceeded as an infinite number of flow lines parallel both to each 


other and to their confining banks. 
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Lelyavskii's work on the Dnieper basin (1893-1904) led, from 
observations Of the velocity and direction of flow at different points 
in the river, to the conclusion that two distinct types of flow co- 
existed within channels. "Converging flow'' causedscour while "diverging 
flow" left channel deposits. Lelyavskii later proposed that a natural 
channel has its adjacent flow lines parallel neither to each other nor 
to the banks, except in the flow zone near the banks. (Kondrat'ev, 
1962) 

Fargue (1908) conducted fluvial studies based on the Garonne 
River. His basic observation was that the maxima of channel depth and 
channel curvature are not coincident, but rather, that the point of 
greatest depth in a channel is displaced downstream from the cross- 
section of greatest curvature. This phenomenmnhas since been confirmed 
(Chang, 1969; Leopold and Wolman, 1960). Fargue also observed analogous 
lags between points of inflection of channel curvature and shallows of 
minimum depth. He concluded that variations in channel curvature 
determine the stability and roughness of the channel such that curvature 
changes would directly induce changes in the channel profile. Greater 
curvatures were expected to produce deeper scour farther downstream. 
(Kondrat'ev, 1962) 

At the turn of the century, then, fluvial theory was already 
exploring the existence of connections between channel curvature and 
depth, and the notion of transverse circulation was available as a 


possible explanation of the initiation of meandering. 


B. Process Analysis 


Another school of fluviologists was contemporaneously concerned 


with the association of river meandering and the rotation of the earth. 
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Eakin (1910) wrote of "unbalanced lateral erosion" due to the earth's 
rotation, and even as late as 1943, Quraishy (1943) championed this 
theory. The theory essentially considered the Coriolis force to be the 
force causing rivers to erode on one bank preferentially, thus initiating 
instabilities which could set up and perpetuate meandering. The effect has 
since been considered, in conjunction with the velocity profile of a 
river, to explain the genesis of secondary flow in rivers, and thus the 
genesis and continuation of meanders, (Scheidegger, 1970; EZinstein, 1926). 
However, the effect has also been given less credit for the cause and 
maintenance of meanders by Chang (1969) who discounted the influence 

of the earth's rotation on the basis of a comparison of the cross-spectral 
analyses of depth and curvature data, using both "relative'' (the positive 
and negative sign convention for curvature of flow) and "absolute" (all 
deviations assigned their positive or absolute values) curvature. Noting 
that the use of absolute curvature gave the higher correlations with 

Other parameters, Chang assumed that the sign of curvature had no physical 
significance with respect to width, depth, or cross-section. 

As the theory of transverse circulation, or secondary flow, 
developed, new schools of thought on meandering grew and gained acceptance. 
Matthes (1941) and Friedkin (1945) typified one school with the consi- 
deration that irregularities in bank erosion and sediment deposition 
constitute sufficient explanation for the meandering phenomenon. It 
was proposed that rivers flowing through erodable materials systematically 
transfer bank material to downstream reaches while eroding upstream 
banks, and the result of the erosion-deposition process is the gradual 


downstream progression of the entire river pattern. Valley-slope, bed- 
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load, discharge, rate of erosion of banks, bed-resistance, and 
transverse oscillations (not secondary flow, but seiche or gravity 

waves) were presented as the hydraulic variables interacting intricately 
to produce the migration. The actual physical interactions were not 
described mathematically, however, and the influential role of secondary 
flow was ignored, (Friedkin, 1945), or rejected emphatically (Matthes, 
1941) by the erosion-deposition proponents. 

It has since been emphasized that even in the absence of bank 
materials, streamflow often easily assumes a sinuous path. Supra-glacial 
streams that melt rather than erode and deposit their bank materials 
(Knighton, 1970) and even streamflow generated and maintained simply by 
water density differences, best exemplified by the Gulf Stream (Leopold, 
Wolman, and Miller, 1964) have exhibited the classical sinuosities that 
rivers flowing through erodable channels seem to imitate. 

Seiche effects were explored further by Hjulstrom (1949). River 
turbulence was introduced as a likely hydraulic component of meander 
genesis, but little turbulence theory was available to provide adequate 
explanation. The seiche effect was given as a component of turbulence, 
while secondary flow escaped mention. 

The use of flumes and hydraulic laboratories led Mockmore 
(1944), Shukry (1950), and Prus-Chacinski (1954) to revive and propound 
the notion of the influence of secondary currents in streamflow on 
river morphology. Spiral motion was observed in straight flumes as well 
as in curved flumes, and was connected with the observation of movement 
of bedload both downstream and laterally, towards the inside of bends. 


Helical flow satisfactorily provided theory supporting a differentiation 


in the magnitude of longitudinal and transverse tangential stresses on 
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che river bed and thus theoretically accounted for both the scouring 
and depositional actions of the river. 

Still, all these classical explanations of sinuosity, proposing 
that plan form meandering is induced by the earth's rotation, regularities 
in erosion and deposition, seiche effects, and secondary currents, failed 
to allow the prediction of river geometry. Attempts to calculate 
meander forms on the basis of these models have met with minimal success. 
The processes definitely participate, in varying degrees, in the main- 


tenance of the sinuous habit of a river, but none sufficiently defines 


the initiation or persistence of stream meandering. 


C. Form Analysis 


Another approach towards the understanding of the mechanisms 
of meandering entails the examination of river geometry and the building 
of models on assumptions of form. Initial attempts assumed meander 
loops to approximate segments of circles. Chatley (1940), Mockmore 
(1944), Shukry (1950), and even Prus-Chacinski (1954) postulated that 
streamflow followed paths constructed of combinations of circle segments. 
Ippen and Drinker (1962) also used arc combinations, but admitted the 
POSsiIDELity Of short, straleny reaches, as well, in order to pereer 
duplicate MoE tM work. Yen (1967) considered a meandering channel 
with circular bends connected by relatively short tangents in order to 
explain spiral flow in natural rivers. 

A more successful model for planform of streamflow is the 
sinusoid, postulated by Leopold et al (1960), Kondrat'ev (1962), and 
Toebes and Sooky (1967). Langbein and Leopold (1966) further proposed 


that the thalweg path approximates a sine-generated curve such that 
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river curvature is a harmonic function of a distance along the thalweg. 
The model is superior since it allows a random element to enter the 
determination of form. The development of the model assumes the 
curvature elements to be normally distributed and to proceed in the 
fashion of von Schelling's random walk (1951). The sinusoid patois 

the model of the river by minimizing angular deviations in the channel 
curvature. At the same time, horizons for analysis of form are expanded 
through this model. Fourier analysis can be applied to the sinusoidal 
geometry to obtain form generalizations not previously possible with 
other models. 

The analysis of visible planform properties of natural rivers 
constitutes the basis of an alternate approach taken towards the develop- 
ment of theory in Piet geomorphology. Bates (1939), Inglis (1949), 
and Leopold and Wolman (1957, 1960) studied relationships among meander 
wavelength, meander width, and channel width. From these and similar 
studies, power functions using these variables have been formulated. 
Quantification of planform was more empirical and at the same time no 
less subjective than ever before in this approach. '"'Typical" sections 
of river were selected for study and definitions of the wavelength then 
depended upon the arbitrary choice of end points of the typical bend, 
(Langbein and Leopold, 1966). Even so, the separate analyses of 
process and form were useful in that they provided the natural and 


necessary background for the inspection of process-form interactions. 


D. Process-Form Interplay 


Jefferson (1902) presented an initial attempt to relate meander 


dimensions to hydraulic variables. He measured mb, the maximum width 


a 
- barbs e ; Ht rn 
oe 


i. =f 
ods 4oaas oF suomi an Ober or avert i? itll a ai laboe 


P in ae 7 | 


«fows att ofs stools ow wank alia His ee 


od? vamvere Isbom dno sergio tovab edt -mugt io ae 


_ 
s43 ntl bsssorq a7 base ae viewed 64 09 ainemels Br 


artioome btoeonte aff . (1201) Slew moboos =? antl fodoe nov to —_ 


lanneds O44 wl enortel vob <thlugns onixtmicin yd travis ort3, to tat boat 6a 


bibmagxo oth cro 15 eleyisns 03 aviee? yon ,siata omee sda JA oma avms > 


inbsosbata add oF botiggs vd iho Bhaylenp aeizuot .isbhg ssp 


oT i4 
oes 
a 
_ 
a 


vistobrs1q don enatiestinrensy mot nisido o3 yrTismesy 


dgitw aidiazog 

i] 
+8 Lebont rod: ° 
atoviy Intuten Yo asisiaqosq mrdtitelq aldtety to aiexiane ed? ie 


-qgolevab wit e251 49 44nd dasosnga sdenvetle we to eiaad sis asatel aes o>. 
rom 


CONOL) aflunT (REL) essed eyeolorgtanolg feivul3 ml goed? Yo 8 
“TeHAEDM, WOME eqisenclis lox herbuith (o7e! SPed) nomlow bas stogood bn 


ne iinla bos oped mort ilo bbw! Penesds- bree fibiw sebracm .figaed 
a 
ad sGhd agoldsinay 5es 1d snteu BEDAS sia Tawed ant ous 


twb [Soielgms s1om anw mreluelq lo soldsak ee ans) 


arolsose: Isotyy?t" SAR atid mi, sieted eve neds avi sok sant 


worntt déanelavaw ody Yo enaiginiteb bee ybeseg ot ba}aa tes Saw 219 


14d. instaqvs.ofa 19 2 ns bus. to apled> yasvatdae ors rio bs be 


70 eoevisnas: Salts Ba ont <4 naval ABR bioqoat bas. rT ~ 
-_ 


ics, bored am od kabivorg vous tools ait teen oo ra bee aae os 


- - 
oe seHOtisex9 Il midt-gasoorg 19. righ angent 9 i smug oe 1822 999% 


ee hoe 


a — um S39bdb3. 0 


-_ a 
ont : Is 


7 


ef tne “meander Delt’ tor a mimber of rivers, and "by listing the 
corresponding "md", or discharge, and the corresponding "w'' or channel 
width, he concluded that the width of the meander belt is a function of 
stream width, and the lateral and vertical cutting of meanders is most 
intense at flood stages. Although his study was subjective both in 

its choice of reaches, and in the subsequent measurement of meander 
belts, his conceptual contribution to the study of the interaction of 
form and process in meandering was important as well as innovative. 

Subsequent studies of geometric-hydraulic relations developed 
and improved upon Jefferson's ideas and methods. Friedkin (1945) 
explored the physical influences of slope on various aspects of the mean-~ 
dering form. He found that increased slope caused increases in meander 
width and wavelength. 

Inglis (1949) rigorously extended the physical characterization 
of streamflow and found that meander dimensions generally seemed depen- 
dent upon the square root of the discharge. Leopold, Wolman, and 
Miller (1964) also accepted the square-root-of-discharge relationships 
as fact. 

Leopold and Wolman (1957) studied the interaction of variables 
in laboratory rivers and concluded that the strongest relationships 
occurred between width and discharge. Discharge was demonstrated to 
be clearly influential in the modification of channel width, and width 
subsequently tended to control changes in meander wavelength. 

Garlston (1965) studied natural rivers to investigate their 
hydraulic relations. He admitted that a major shortcoming of much of 
the fluvial work conducted (including the previous work of others as 


well as his own) was the subjectivity inescapable in the choice of 
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reaches examined. Often, workers based their findings on the measurement 
of a eines meander rather than on a statistically valid sampling of 
meanders, and choice of workspace depended too highly on the availability 
of appropriate data, which was scarce. Carlston contributed a new concept 
to the theory of meandering: the dominant controlling discharges in 

the modification of meander wavelength are those ranging between the 
highest monthly discharge for a given year and the mean annual discharge. 
This conclusion was offered to replace Inglis' notion that discharges 
slightly greater than bankfull flow were the dominant modifiers of 

river meander wavelength (Carlston, 1965). 

The proponents of this form-process interplay approach to the 
investigation of meandering generally agreed that linear planform 
parameters such as wavelength meander belt width, channel width, and 
wetter perimeter, depended directly upon the square root of the river's 
dominant discharge. Two distinct definitions of dominant discharge were 
involved (Inglis, 1949; Carlston, 1965). Slope was also accepted as 
influential on the basis of Friedkin's 1945 lab studies at Vicksburg. 

The proponents succeeded in their attempts to bring rigor into the 
discipline. However, fluviologists generally examined pianform on the 
basis of downvalley distance and width of meander belt, ignoring actual 
stream path-length. Meander belt was measured as the transverse 
distance between tangents to two adjacent or appropriate but opposite 
river bends, and it was used to describe the magnitude of the stream's 
oscillatory activities. Wavelength was approximated as the downvalley 
distance characteristic of the river's meandering wave. The path-length 


was not yet given any significant theoretical importance. 
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Mackin (1956) and Schumm (1963, 1967) were among the first to 
examine the role of channel pattern in the form-process interaction. 
Mackin (1956) found the relation between sinuosity and cross-section 
to be such that as sinuosity decreases from tortuous to straight 
channels, the width-depth ratio of the channel increases. Sinuosity 
was then defined as the ratio of the path-length of the stream to the 
down-valley distance, a parameter still necessitating the inclusion 
of subjectivity into theory, because for all but perfectly regular 
meanders, the end points of the down-valley line must be selected 
arbitrarily. 

Schumm (1963, 1967) compared the same sinuosity parameter to 
parameters of bank cohesion for several American and Australian rivers. 
He found that greater sinuosity values tended to characterize streams 
with more cohesive bank materials. Schumm defined cohesiveness by the 
percentage of silt and clay found within the bank material. 

Considerations of stream sinuosity opened a new theoretical 
dimension in fluvial geomorphology. Different river patterns could be 
investigated to help strengthen previous notions and create new theory, 
but the elimination of subjectivity in the parametric characterization 
of sinuosity remained a problem of paramount importance for fluviolo- 


gists. 


Ey. COntributions of Leopold: The Sine Generated Curve 


Leopold and Wolman (1957) introduced a new and substantially 
superior definition of sinuosity when they considered it the ratio of 
the thalweg path distance to the valley length. As the thalweg is the 


river's own "preferred path", in that the thalweg will actually be the 
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river's course at very low discharges, the new sinuosity parameter 
described the river's meandering activities much more accurately. This 
concept was closely related to another significant discovery: the 
thalweg of a straigiit channel exhibits a sinuous tendency as it flows 
between pools and through riffles. With some laboratory experimentation, 
Leopold and Wolman (1957) demonstrated that a straight channel did not 
necessarily, or in fact usually, possess either a uniform stream bed 

or a straight thalweg. 

The model thalweg inferred from the laboratory work had a wave- 
length of L x 2/7 radians (one complete sine wave) and was characterized by 
relatively deep pools at bends with shallows at points of inflection. 
This was offered, then, as a finding common to streams of both sinuous 
and straight channels. The wavelength was generally a function of 
width and so indirectly related to discharge. Bankfull discharge was 
assumed the best flow parameter in deriving flow-planform relationships. 
(Leopold and Wolman, 1957) 

Leopold, Wolman, Miller, and Brush (1960) conducted laboratory 
experiments that produced beautifully sinuous channels containing no 
straight sections. They explained the meandering tendency as an energy 
dissipation technique. If bed material was of insufficient size to cause 
bed stabilization, the river would be free to create random bank 
projections. Local chance deflections would set up systems of eddies 
and vortices which would strain and erode the banks in the streams 
process of effecting internal energy losses. 

Both sinuosity and energy considerations were greatly expanded 
in Leopoldand Langbein's 1966 classic work, "River Meanders''. Random, 


or stochastic, processes were finally allowed in meander development 


pes net uae 
Ree ipeae Mea 
mm. bee : issn ant RN 
ena shit tad 
sen neh Aen 


ros Daria 


ats Wit: Tian sehgtaean dy bi ne 
bad aonze, penny ‘tanto meadion a 2 bk 
sovow é bai drow wrodetods? add ei elec eee ar 
4d bogtaqaoo dla sew bag Covew ote sdalgnos sap) andi ber we x ae 7 ty 7 


waltoodtri Yq soetdn mm, ewoP fade sist obnsd tn efocg 

anamia tized To nets of Yom gaehaid 4 26 .neds bosso gawe 

te tplso0 « yf lerney anew iigieliver st? . dahaiaita, naa 
any oytivieth Dutfeed .. sanrivlontb ot) booates aobr bal oe Om eee 


er, 


aintpeetseths. nilislaswoltgalvirsh ot cotenwuig’ wl? deed sa: | 
| (SMP aoe ious : 

weeihseded bstoubaon (G0!) fous le +,cal lan ie ow bleqoxt | a 

Oo sotalam alguns’ auopinta iigtetead IwonGosy dads Sram Teg Ke 

YRiAs As es vonsbas? ainhwating ce a toclalgne wT enol seay saghaneta 

mine OF 9ate deoiolyivant Op jew lelesdee hol TY sabiphovlzey nab seqeaanb: 

dned, mobyna. si941.03) gost of ‘Hleow yoyni4d erly ee ' 

valde Ye dashes! es tow bivow enofaopltab sanptts Iepat eno) somo 


dansyi. avid «tt ashied cdi ebesy tue Witte diuew stoide 


ssnan Pe stgeliel dames galtectee Se 


| | beers sigty. ns ee deaneraie Retdaad An vane de as. 
: _ a es he oe 
camnaiiints lint {beg 


ar a | , 
: ” 7 oe 


16 


theory. The random model used was based on von Schelling's (1951) 
random walk in which successive moves on a surface were independent 
with respect to direction travelled each move, but were governed by a 
probability function such that, relative to downvalley direction, the 
expected directional deviation at a point was a sine function of the 
channel distance travelled. 

Streamflow was likened to the chance-directed particle which 
should supposedly seek a path eventually minimizing variation in direc- 
tional changes. The sine-generated curve was postulated to be the flow 
pattern best allowing such minimum variance (given symmetrically curving 
paths of equal length), thus minimizing total erosion, as it was shown 
that banks are eroded at a rate directly dependent upon channel curvature. 
The sine-generated meandering was assumed, on the basis of field work, to 
provide the most uniform water surface, signifying maximum uniformity 
in energy expenditure per unit of distance along the channel. Thus, the 
sine-generated curve model was shown to eliminate concentrations of 
energy expenditure, following Newton's second Law of Thermodynamics. 
(Leopold and Langbein, 1962; Scheidegger, 1967) 

Leopold still considered wavelength as a downvalley distance, 
however, and thus retained the subjectivity of previous analyses of 
form. A method that would eliminate this subjectivity, recurrent in 


hydraulic analyses of stream activity, had yet to be found. 


G. Speight's Meander Spectra 


MacKay (1963) initially proposed the application of the spectrum 
to studies of fluvial meandering, but he did not apply the technique 


himself. Speight's work initiated the actual application of the 
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technique. 

Speight (1965) sought to overcome three significant biases in 
fluvial research. He complained first of the subjectivity involved 
as workers chose "typical" meanders on which to base wavelength 
statistics, and then of the error inherent in measuring wavelength 
along straight lines that probably have little relevance to actual 
flow direction. He finally concerned his research with the tendency 
of previous authors to ignore the possibilities of multiple wave- 
lengths in stream form. 

Speight used sets of air photographs of the Angabunga River in 
Papua in his research on curvature. He calculated successive angular 
deviations at points measured along the thalweg, separated by a constant 
interval along the path. He then applied spectral analysis, a technique 
described in detail in Chapter 2, to the derived series. The technique 
is used to separate complicated sinuous series into Fourier components 
assigning each Fourier term its variance contribution within the whole 
meandering series. Speight used the variance statistics as objective 
parameters of sinuosity. Multiple wavelengths were described for the 
first time through the use of this technique. 

Attempts to relate the significant spectral peaks, or wave- 
lengths, to tree-to-tree and grass-to-grass river widths were largely 
unsuccessful. The shape of the derived spectrum seemed to suggest the 
influence of channel cross-section, though, and this observation again 
opened new dimensions in fluvial work. Speight then attempted to 
compare spectral analyses on curvatures derived from the same river 
but from different years of flow activity. The air photos provided 


24 years of the necessary data. No orderly meander migration was observed 
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but meander wavelength seemed constant. Spectral analyses was shown to 
have significant potential for the analyses of many other hydraulic 
variables, and Speight encouraged readers to expand research with this 


too kand(Spereht, (9650. 


H. Later Uses of Spectral Analysis 


Thakur and Scheidegger (1970) continued the use of stochastic 
models in the study of sinuous river behavior. They considered a river 
course to be a stochastic line subject to certain geologic conditions. 
Spectral analysis was used to examine the curvature series of the 417 
chain models they generated in an attempt to discover an expected or 
most probable spectral value. Although the expected value notion later 
proved less significant than they had hoped, Thakur and Scheidegger 
were able to conclude that the frequency distributions of deviation 
angles was reasonably Gaussian, an assumption of spectral analysis 
which had been previously unproved. They also developed the use of 
transition probability matrices for the determination of the mode or 
pattern of river bends, which was reflected in changes of sign in the 
series of the deviation angles. 

Scheidegger (1970) in Theoretical Geomorphology, cites spectral 
analyses as an important tool in the measure of "wiggliness" of natural 
lines. The substantial potential for spectral analysis in turbulence 
studies is also recognized in his work. 

Spectral analyses of sand waves has been very successful and 
shows great potential for the quantification of descriptions of bed 
roughness. Nordin and Algert (1966) were able to distinguish the 


spectra of sand dunes within long records, and they proposed a second- 
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order Markov process for wave formation. Ashida and Tanaka (1967) 
further developed theory on the second-order Markov model, and they 
used spectral analyses in the determination of the propagation velocity 
of sand waves. O'Loughlin and Squarer (1967) and Squarer (1968) 
recommended the use of spectral analysis to standardize descriptions of 
sand and waves. Nordin (1971) conducted flume studies of sand waves 
and concluded that the spectra of the waves were similar regardless of 
scale of phenomena. The general shapes and movements of both dunes and 
ripples in large ne small flumes showed similar characteristic profiles, 
exhibiting only scale effects, when present, to provide differentiation. 
Nordin's findings agree reasonably well with Yalin's work of the same 
year. (Nordin, 1971; Yalin, 1971) 

Chang and Toebes (1970) extended the use of spectral analysis of 
river curvature to the examination of absolute curvature of meanders. 
For this they used only absolute values of angular deviations along the 
Wabash and White Rivers in Indiana. They recommended the use of 
absolute curvature in cross-spectral analysis, although objective reasons 
for this were not clear. Chang found that the use of absolute curvature 
simply allowed higher correlations with the other parameters studied, 
and for this reason, he chose it as the superior statistic. Meander 
planform was demonstrated to be controlled by geology and soil, and 
flow rate and channel slope, through the use of contour maps of 
spectral amplitude statistics for overlapping successive downstream 
reaches. 

This method, the "moving spectrum", is a way of avoiding serious 
non-stationarity in their data. Chang and Toebes were able to diffe- 


rentiate glaciated from non-glaciated fluvial landscapes on the basis 
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of their spectral findings. Meander wavelengths in unglaciated bedrock 
tended to be longer than those of glaciated regions. 

Chang (1969) extended previous work, using the White and Wabash 
Rivers, to include cross-spectral examinations of width-depth, depth- 
cross-sectional area, width-cross-sectional area, curvature-width, 
curvature-depth, and curvature-cross-sectional area. Much of his work 
corroborated earlier theories, such as those relating planform 
parameters to the square-root of discharge, and others observing an 
increase in wavelength in the downstream direction. The analyses of 
Chang (1969) and Chang and Toebes (1970) constitute a long anticipated 
advance in objectivity and rigor of analysis not attained in earlier 
studies conducted by other fluviologists. However, a variety of channel 
types must yet be examined, and comparative studies must still be 
conducted in order to test the findings of Chang's work based on the 


Wabash and White Rivers. 


I. Other Recent Theoretical Contributions 

Church (1972) applied various spectral analyses to sandur 
rivers on Baffin Island. He found a dominant pattern of peaks to be 
consistent, recurring at short, intermediate, and long wavelengths. 
Church found dominant wavelengths, expressed in terms of channel width, 
to persist as multiples of es 

Yalin (1970) used theories of river turbulence to explain the 
occurrence of sand waves and the regularities of meanders. Yalin's work 
predicted dominant meander wavelengths, expressed as multiples of mean 
channel width, to occur as integer multiples of 2%, and dune wavelength 


to be similarly proportional to 2mh, where h is the mean water depth. 
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Further, Yalin has made a most important observation: that, 

a systematic disturbance caused in the behavior of the 

largest macroturbulent eddies by a permanent discontinuity 

must be accompanied by the same and the reverse kind of 

disturbances occurring alternately, with equal length 
intervals, along the direction of flow. The dunes and 
meanders are the developed versions of the alternate 

erosions and accretions caused by these disturbances in 

the velocity field, and thus in the transport capacity 

of the flow. 

(Yahin, -1970) 

Surkan and VanKan (1969) extended work on random walk models, 
incorporating hydraulic considerations. They found existing models to 
be insufficient, in that simulated paths lacked the regularity of 
oscillation and freedom of erratic behavior found in nature. The 


inclusion of a gradient bias in random walk models improved models of 


meander patterns substantially. 


J. Conclusions 

Spectral analyses of both planform and hydraulic characteristics 
of meandering rivers provide objectivity and rigor that previous 
analyses lacked. The stochastic approach to river behavior provides 
a superior form of analysis at this stage in the theoretical development 
of the meandering intrique, Stochastic models facilitate an understanding 


of process and form in fluvial geomorphology. 
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CHAPTER II 
STOCHASTIC PROCESSES AND THE VARIANCE SPECTRUM 


The major classical and contemporary explanations of the initiation 
and persistence of fluvial meandering were briefly explored in the previous 
chapter. Among the processes considered to be meander-generating were 
activities of secondary currents in the cross-sectional plane of a stream, 
influences of the earth's rotation, regularities in bank erosion and 
subsequent deposition of materials, and seiche effects analogous to lake 
seiches. Attempts to calculate and predict meander forms and behavior 
on the basis of these models have met with minimal success. The processes 
examined definitely contribute to the establishment of the meandering 
habit of rivers, but none sufficiently explains either the genesis or 


geometry of meandering. 


A. STOCHASTIC MODELS 

Davis! (1893) firmly established mode of landscape analysis 
emphasized development according to a well-determined progression of 
stages. This leading role assigned to historical context in the 
examination of developing landforms was conceptually significant in that 
it allowed sensible development of theory on a macro~scale. Valuable 
generalizations were produced and some processes were isolated for 
further consideration. But studies made on a more micro-scale, 
examining detail within such a deterministic framework of logic, were 


complicated by the characteristic interaction of a large number of 
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inadequately understood variables. Fluvial behavior proved too dynamic 
for such "rigid" analysis of stages. Details of interactions were not 
satisfactorily understood using this "closed-system" approach. 

The notion of landform equilibrium was inherent in the Davisian 
Stage analysis of landscape development. Other disciplines, such as 
biology and physics, commonly examined equilibrium interactions of 
ensemble variables in terms of systems theory. Chorley (1962) discussed 
the value of the application of systems theory to work in fluvial geo- 
morphology. While the work of Davis and others in fluvial geomorphology 
was influenced somewhat by closed system logic, Strahler(1950) assumed 
that natural phenomena operated in both closed and open system fashions. 
Chorley (1962) recognized Strahler's work as a great theoretical advance 
for the discipline and he encouraged the development of theory on the 
assumption that natural systems are open. 

In systems analysis, a closed system is an assemblage of objects 
and their relationships. Boundaries across which no energy export or 
import occurs are required in the definition of a closed system. Changes 
in the system are deterministic, progressing towards the establishment 
of equilibrium conditions. The progression towards equilibrium in a 
closed system simultaneously entails the maximization of entropy within 
the system. Thus, potential for change gradually decreases towards 
system stabilization. In this way, assumptions provide that initial 
conditions of the system are sufficient to determine the eventual 
equilibrium situation. While the deterministic model contains much 
useful information, it lacks the flexibility attained when systems 
are considered to be open. 


Open systems assume assemblages to be characterized by dynamic 
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energy inputs and outputs, thus including closed systems as the special 
case of zero energy exchange. The theory of the open system's steady 
state allows more flexibility than the closed system's idea of equilibrium. 
Dynamic equilibrium, the steady state of open systems, is the dynamic 
condition in which the open system continually adjusts to externally 
induced physical changes by changing its own form. The concept of adjust- 
ment allows a more detailed analysis of form-process relationships than 
previously possible, and it admits a more liberal view of changes. The 
Open system more easily directs attention to an entire landscape assemblage 
as it acknowledges the complex multivariate structure of most natural 
phenomena. It is thus a more realistic model for the study of river systems. 
The most obvious advantage of the open system approach is its 
scope for the consideration of natural new inputs and their effects on 
the system. Random processes are thus allowed as stimuli, and responses 
to these are considered important. The model that best explains the 
behavior of hydrologic systems must recognize that the processes involved 
are not entirely deterministic but can be characterized as stochastic. 
Precise descriptions of river morphology and behavior are valid 
for a river only at some specified time. The dynamic nature of fluvial 
processes is manifest over time as extensive changes are induced both 
gradually and catastrophically. Such changes become reflected in both 
the river's planform and bedform through adjustments by lateral planation, 
downcutting, or aggradation of the stream. The changes imposed are 
predominantly determined by the energy contained in the stream, 
Fluvial processes must be considered as functions of both time 
and space. The analysis of river behavior with respect to such factors 


as discharge, slope, grain size of bed material, sediment load, or cross- 
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sectional form cannot be conducted independently of time and space 
considerations. Since flow constitutes a progression in both time and 
space, the factors studied will vary in such a way that they are inde- 
pendent of distant past or future changes, but dependent upon immediate 
past events, utilizing a portion of the information contained in an 
event of the immediate past. Connections with other points in time and 
space thus contain an element of determinism, but they also allow 
probabilistic relationships. Stochastic models provide for the co~-oper- 
ation of both deterministic and probabilistic processes, and they are 
thus the most suitable models to describe river behavior. 

Figure 2.01 is useful in the clarification of the concept of 
astochastic process. Suppose a man begins walking from a point A, intent 
upon reaching some point B six steps lower and twenty paces ahead. If 
he is permitted to move only forward or downward, then he has two choices 
each time he takes a step. Since he must reach point B, his decisions 
are influenced by his position relative to B, and probabilities may 
be assigned his choices. He has seven possible downward moves from which 
torchnouse (091° 2, 2, 4, 5, or 6 steps down) ¥ and so for any positian 
(x,y), his probability of choosing a downward step is y/7. Similarly, 
he chooses a horizontal step from position (x,y) with a probability of 
by fan Ie Probabilities of owes decrease as the man approaches B (or, as the 
x and y values decrease). The probabilities vary in a manner only 
partially determined by an determining last and next moves, and probabi- 
listically independent of other past or future events. An element of 
influence can be associated only with time-~adjacent moves. This type of 
model exemplifies the random walk and may be represented as 
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where Sr represents the position of the man taking a random step Xi 
at time t, so that Xi indicates the distance covered in any one step. 


N 
Here > Xi may be used to describe the path distance between points A and 


ita! 
Be 

In the sequence of movements, the variable Xi at time Ce 
independent of the entire previous set of X's. The future cannot be 
completely predicted on the basis of the past record, but it relies both 
on the existing record and the probability of the next move. The 
probabilistic element precludes the use of deterministic models in 
recreating or simulating a series produced by the random walk model. 

Von Schelling (1951, 1964) developed a variety of random walk 
models that were later applied to studies in fluvial geomorphology. 
Leopold and Langbein (1962, 1966) used the Von Schelling elliptical 
path model to explain their observation of the tendency of river 
meanders to assume geometrics minimizing the variance of the angular 
deviations along a path. Scheidegger (1967) further developed the model 
to extend his thermodynamic analogy between the movement of gas molecules 
and river behavior. 

Stochastic processes may be used to explain the random nature of 
phenomena changing with time. The stochastic series may be represented 
as the family of random variables X(t) where t, the time parameter, is 
an element of the set of intervals, T, or as: Pesene té T| » the 
stochastic model accounts for both the probability and sequential relation- 
ship of the random variable. Conventional deterministic methods require 
the assumption that the process does not change with time or spatial 
variations, which is inherently contrary to hydrologic observations. 


Any hydrologic time series may thus be represented as 
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6c 7 At, where xX. is the Series, 2 Lene 
deterministic component, andé t is the random or probabilistic element 
(Matalas, 1967). The uncertainty generally observed in the series is 
thus accounted for byes. and this residual element is assumed always 
to be present. The deterministic element may not necessarily exist in 
the series, and when present, its importance in the prediction of xX, 
may vary from predominant to negligible in the determination of the 
process. The two elements may be mistaken for each other when records 
are very short, and it is essential that time series analysis minimize 
this possibility. 

The amount of information retained or repeated by subsequent 
points in the time series may vary. Dependence between adjacent 
observations may be measured by the autocorrelation coefficient of the 
series. Markov chains are a special class of stochastic processes in 
which the value of X, at any set of times t, depend on the pale’ x, 


hs 


at any immediately previous set of times Le The outcome of any event 

is considered to depend on the outcomes of the immediately preceding 
events, but is independent of all previous events. A Markov chain 

"of order n'' specifies that in a sequence of events, the outcome of each 
event depends on the Bin fons of the n directly preceding event(s), but 
only on these. First order Markov processes are of particular hydrologic 
interest, as their patterns are recurrently found in river series. 
Adamowskii and Smith (1972) found that the stochastic component of rainfall 
generation was adequately described by a first order Markov process. 
Markov processes without autocorrelation of series elements are described 
by the random walk model. 


Any stochastic model by its very nature is not supposed to give 
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a perfect fit to the record it models. The actual record obtained is 
assumed to be only a sample of the many possible samples that could 
occur according to the probability distributions proscribed. An adequate 
model is rather one which preserves the desired parameters essential to 
the interpretation of the data. Markovian models preserve the lower 
order moments of a time-sequential distribution (mean, variance, and often 
skewness), and the first or higher order autocorrelation coefficients 
(Matalas, 1967; Rodriguez-Iturbe, Dawdy, and Garcia, 1971). The preser- 
vation of these parameters is generally adequate for work in fluvial 
geomorphology. 

Stochastic processes are referred to as stationary when the 
statistics of the series are not affected by a shift in the absolute 
time or space origin, and so are reasonable approximations of the sta- 
tistics of the appropriate ensemble. For instance, if the series is 
shifted K time units with respect to the origin, then, for the parameter 
xy 
f(z; t) = f(x; t + k), must be true for all terms, K, under the assumption 
of strict stationarity of the series. A series is called weakly stationary, 
or stationary in the wide sense, if the expected values of the first 
and second moments remain constant when the time origin is shifted, and 


the autocorrelation (Ry) depends only on the time differences (t,-t,) as: 


E [Xe] ae = 4g constant 


and Cov Rear XY ok 


Stationarity is an assumption theoretically fundamental to higher order 
tests made on the time series, but it is generally assumed only in the 
wide sense, when considered at all. 


In practise, rigorous tests of data stationarity are seldom run, 
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Ob ig tum, seldomadesdribeds | Tess for 'stationarity GE %data are “nor 
included in any of the literature to date applying spectral or cross- 
spectral analyses to the problems of fluvial meandering or bedform 
configuration. Tukey (1961) approached the problem more optimistically, 
encouraging a less-inhibited, while cautious, use of the spectrum: 


» . - I have yet to meet anyone experienced in the 
analysis of time series data . . . who is over- 
concerned with stationarity. All of us give some 
thought to both possible and likely deviations from 
stationarity in planning how to collect or work up 
data, but no one of us will allow the possibility 

of non-stationarity to keep us from making estimates 

of an average spectrum .. . The fact that the 

spectrum is changing with time ... need. not make 

it unwise to estimate one, or several, average spectra 
- « - « Once we admit that we are estimating an average 
spectrum, we have admitted that there may well be other 
relevant characteristics of the situation beyond the 
spectrum, that estimation is not complete specification. 
Such an admission . . . is a good thing rather than a 
bad one. (Tukey, 1961) 


It appears, then, that while stationarity should be considered in time 
series analysis, its absence should not alone prevent the use of spectral 
analytic techniques. Degree and form of stationarity of data should, 
however, be remembered in the interpretation of results. 

Three main models of time series are used to simulate stochastic 
hydrologic processes. When the nonrandom portion of the series can be 
attributed to storage in the hydrologic system, simulation is best 
accomplished by the autoregression model, which may be represented as 


SES 


x. = E(x paQ? +e: Xtep : 


t ve 


where x, is the process under study and K is the time shift parameter. 
(based on Chow and Kareliotis, 1970) A special case of the autoregression 
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model is the linear autoregression model of the nth order: 
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mere ays 5 2 6 « are the serial correlation coefficients of the Variat 
x. When n= 1, the process defined is the first order Markov process: 
Se ex py ap C,- Figure 2.02 presents the correlogram for the first order 
Markov process, and Figure 2.03 shows the spectrum of this series. 
Another special case of this model is the second order linear auto- 
regressive scheme, in which two correlation terms arise: 

we 84 Xe-7 7 89 7 * St 
The scheme is similar to the first order process, but shows a slightly 
longer memory, such that events in the present rely on two events in the 
immediate past than on the preceding event alone. 

A model useful in the representation of carry-over in a hydrologic 
system due to water retardation characteristics of the watershed is 
the moving average model. The model may be defined as: 

xX, = ay G. = ay Gen ueusee: as, rae 
where aj, ays + + + a, are the weights and m is the extent of the moving 
average. By virtue of the moving average on the 6's, the model is not 


purely random, but is stochastic. (Chow and Kareliotis, 1970) 


Oscillating variations in the sequence may be represented by 


N a he 2 ee 
the sum of harmonics model, where: x. = 2 (AieosOOT I eketsit FT ) 
; ct 2 aye 
where A; and B; are the amplitudes of oscillation, tr is the period 
Ge cyelteity with | — 1, 2, . + .; and N is the number of record intervals. 


The model assumes the oscillations may be represented as sinusoidal. 


Diurnal, seasonal, and tidal periodicities in hydrologic systems may be 


approximated by the model. The model may also be applied to the meandering 


behavior of streamflow (Chow and Kareliotis, 1970). 
The particular choice of model is not an obvious one, and may be 


aided by the use of the correlogram, as figure 2.02 shows. A correlogram 
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LAG K 


AN INFINITE SERIES OF PERIODIC 


CYCLES 


PERIOD OF RIS NOT NECESSARILY 
THAT OF ORIGINAL SERIES 


AN INFINITE SERIES GENERATED 
BY AN AUTOREGRESSION SCHEME 


DOES NOT VANISH BEYOND ANY 
CERTAIN POINT 


AN INFINITE SERIES GENERATED 
BY A MOVING AVERAGE OF “M” 
TERMS 


CORRELOGRAM TRACE IS A 
DECREASING FUNCTION TO “M” 
LAGS, ZERO BEYOND THE MTH LAG 


CORRELOGRAM TRACE OF 
RANDOM* SERIES, SHOWS 
PERFECT CORRELATION TO ZERO 
LAGS, ZERO CORRELATION FOR 
ALL HIGHER LAGS 


* MEMBERS OF SERIES INDEPENDENT 
OF ONE ANOTHER COMPLETELY 


Re SIGNIFIES THE CORRELATION COEFFICIENT FOR ANY SPECIFIED LAG “K" TIME UNITS 


CORRELOGRAMS oF TYPICAL SORTS OF TIME SERIES 


(AFTER FRIEDMAN) 
FIGURE 2.02 
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C) FIRST ORDER MARKOV PROCESS, 


HIGH FREQUENCY CONTRIBUTIONS 


SPECTRAL POWER 


0 all a2 pis) 4 25 s 


=—— 


SPECTRAL POWER 


0 ail ne a3 -4 HS 


B) NORMAL WHITE NOISE (RANDOM), 
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D) FIRST ORDER MARKOV PROCESS, 
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E) SECOND ORDER MARKOV PROCESS, F) THEORETICAL COHERENCY 


(INTERMEDIATE FREQUENCY 
CONTRIBUTIONS) 


SPECTRUM FOR BIVARIATE 
AUTOREGRESSIVE PROCESS 


SPECTRA OF TYPICAL SORTS oF TIME SERIES 


(AFTER JENKINS AND WATTS) 
Dishing 2.0.3 


diagrammatically shows the effect of a time origin shift on the 
autocorrelation. ie Series 16 correlated with itself using a variety 

of lags, (denoted "m'" for a series of "n'' points) so that a comparison 
of the fit or correlation of the series to itself when shifted temporally 
or spatially may be obtained. The correlation is thus highest when the 
Series are tdentical, as at"a tas "or zero units-of time. The auto- 
correlation increases as the series approach identity. The correlogram 
presents the effect of each possible lag on the autocorrelation. 

Every stochastic model includes 6,, the element of chance or 
probability. In any dynamic system, chance perturbations or fluctuations 
must occur and cannot be predicted. The element of uncertainty makes 
the use of deterministic models unsatisfactory for the study of hydro- 
logic series. The elements of change over time and space in the behavior 
of fluvial systems are far too important to be ignored. The most 
reasonable approximation of the actual processes must be attained through 


the use of the stochastic assumption. 


B. TIME SERIES ANALYSIS 
Recognition of the influence of time and space on series 

describing river behavior has led to the recent introduction of time 
series analysis into the discipline. Time series refers to a set of 
measurements recorded sequentially in time, with measurements usually 
sampled regularly from a continuous record (recorded at a constant time 
interval). Hydrologic series are constantly undergoing changes in time 
or space by natural or anthropogenic processes, but the changes may not 


be obvious in the actual record, or they may be difficult to detect 
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Haan (1972) recognized that one recurring problem in hydrologic 
studies is the lack of adequate data. The obvious solution to this 
problem is a systematic collection of long records of appropriate data. 
Since. long series tend to contain information about the mean, and thus 
give better pictures of the processes under study (Matalas and Langbein, 
1962), such a solution would seem imperative, but other problems are 
raised. If short records can mask the nature of the changes they contain, 
how much more confusing the interpretation of reams of continuous data 
will be! Methods for the analysis of long records of continuous time 
series are essential to a better understanding of the information they 
contain about both form and process. 

The time series analysis model most suitable for the eventual 
explanation of meandering river behavior must be able to detect persis- 
tence in the data in all its forms of occurrence, and it must be sensitive 
to time dependence or serial correlation of the data. Meander and bedform 
data are likely to exhibit some combination of the four major components 
of time series. A periodic component, or one that exhibits regular 
reversals of curvature (approximating a sinusoidal form) is expected to 
characterize the data. A completely random component, or one showing no 
recurrent pattern, will likely exist within the data, as well, obscuring 
other process patterns. Many temporal and spatial series exhibit a long- 
term trend component so that the series changes systematically through 
time. Finally, the sample may encompass a major change in the series 
which is not really periodic; often incorporating some random fluctuations 
in period, phase, or amplitude; but is quasi-periodic or oscillatory in 


nature. Figure 2.04 illustrates these possibilities. 


mietbed bae tehiseM octal oad 16 udite ata Isiv45 26 romalinoged, wait oe: 


~~ te Aa 


o i 
4 f, 
. 


ees 
siyuivabyd wi « dm, leas 
nbiia 03 03 tte ie aia pirat 

-sieb eseimyornge Io ar Ye watszeltan seine 
guilt ben adm wih anaes wiaauas a 
ndodgnet bap colsset) ybuss zabnu, esxqos02q) 2a do —, 
um nmetdor soto Jud oviibreqml mesa bivow notivloa 
~HisIi95 yard eagnnna a9 lo Sssugon odd daem a62) ‘thong 
wiab evauntsnos 26 gmk ST Poon dtonteag vant ~aild-gorietineo wan eae at 
sini} avoutlgans Yo ubsosax gnot Fo pieyidus of2 ot sbodism | 


yea nolsgmaatiit ots Io gribnpievobru waded B03 [attnoges otn 
seednnnen oda tol dod dugde pbeamo . 

feutrovs off 10% sitbotue tegm ishom aleglone eelave emia sat <a 
-sleinq Jodtab 69 aida of saum azolvsled aevts griasbnesm Jo ete a B if 
ovittusee ed taum ti hag .oonstiw2a0 to emto? esf [fe n} @a8h amd, a ran 


#onsfoqnos rwofew toi ads 16 rofatinideco oaoe Didliee os gladil sag enh 
telisgoy euididay ded’ sito vo ,j¢eneq@oo o}borxag A .eshase on lo a 
i : 7 : - 

ed bodoeqxe ci (io) Iebicsunla s oi Jemixeazges) sab3evau5 ee eS 


on gaiwota “sno uo ,faenegqmey wnlinas ¥in {odin A .a20h of ssliioetado: 
“ai 

‘ 7 
-gyol w didides aifsee Untazqe ban Jorqmed YoeM *samegdeq sunconq aild9" : 
Wgitoint yilpolsemeteve saginds aslise ‘odd Jodi 08 1Xeeoyeeo baer 
ow ts eee 

gobi Sdi nt vanes iol sw. wesqhone Yam otquse ans, geet ork I: 
ancltsuioul odbiat singe Ani tusequosnt codIto i yiteds Yor a: 
nt yroasl ioe 1% eibpiteq-tonp et wd ysbosl lqam ae sendg 


solailidieaiq aeoity aoasccbev th dm.e: 
y 


gilavosdo , flow ae , 6306 afd nifjinvet<ixe qisdil Liew .ayasieg 3 


A) PERIODIC 


Ava 


B) COMPLETELY RANDOM 


C) LONG TERM TREND 


D) QUASI-PERIODIC 


COMPONENTS OF TIME SERIES 
FIGURE 2.04 


Each of the four possible components requires special treatment 
at some stage of the series analysis. A long term trend, if present, 
is usually removed at the beginning of the analysis by some sort of 
regression or moving-average process. Some forms of analysis instead 
consider a trend part of a low frequency oscillation, since "long-term" 
is a relative concept. A trend is usually assumed to be independent 
of the other components so that its removal does not affect any of the 
series' moments. 

The general assumption of stationarity, mentioned above in 
section 2.A. and again in chapter 4, necessitates the removal of any 
trend present in the data. Stationarity in the wide sense requires the 
expected values of at least the first two moments of the series not to 
vary with origin displacement, while the autocorrelation function is 
required to depend only on time differences and not on absolute time. 

If higher order moments are similarly dependent only on time differences, 
the series is said to be stationary in the strict sense. If strict 
stationarity holds for all series in the appropriate ensemble, the 

series is further said to be ergodic. Available statistical methods 

are designed to analyze stationary series, often first assuming 
ergodicity. If the series cannot be assumed even to be weakly stationary, 
it should be transformed to a stationary series before rigorous analysis. 
Torelli and Chow (1972) give a detailed discussion of the problem of 
Stationarity in hydrologic time series. 

The most important aspects of time series, when the objective 
of study is discovery of the underlying physical mechanisms, are their 
oscillatory and periodic tendencies. Oscillations are commonly irregular 


in character, exhibiting erratic changes in amplitude, phase, and period 
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of the patterns. The major part of the theoretical work in time series 
analysis is concerned with obtaining an approximation of the variation 
in the series. This was originally approached through examination of 
the correlogram, but confidence limits and significance tests have never 
been successfully applied to the correlogram. For this reason, the 
technique was partially abandoned in favor of the analysis of the actual 
variance of the series. The correlogram is still a useful tool in the 
classification of processes underlying the series, but it is not desirable 
as the final product of the analysis since results cannot be easily 
compared in the absence of confidence limits or significance estimates. 

Periodic analyses commonly use sinusoidal models as shown in 
Figure 2.05, where: 

20 

y =(4 cos Loot ©) 
when A represents the wave amplitude, @ represents the phase angle or 
shift from the origin, and T represents the wavelength. The early 
assumption that time series were simple sums of a given wavelength led 
to harmonic analyses and later to periodogram analysis. The periodogram 
plotted periods present against their squared amplitudes. It was found 
to be a very unstable measure of periodicity of real time series, 
possibly due to random fluctuations in the data, and was then abandoned 
in favor of variance analysis. 

The variation of a series may be calculated as a spectrum of 
contributions to the variance made by a variety of frequencies of oscil- 
latory waves. The contribution to the spectrum by any frequency is called, 


by communication engineers, the power of the frequency in the series. 


The technique used to compute the spectrum of the time series is known 


(0 
LO 


ot yadeset ats 307 miss 
Javton 614 10 elevison qo ee onan mare ig . 
sd3 nt Inox [widen s thaag al mssgoloas6> ad estise sa Se sonsitsv 2 
sldaxtash tom al 3k jud ,2Slasa. ‘ila gay fixabaw Gutman | sellcngeis aaah } 
yitass sd jones? esfrest somber: etey!sne Ri aia Jouborg fant? as oe 
| ,asigmite® Soanolliggre’10 pate} | sities Ao ‘hcligat orfa mt besaqmoa 


| . idee © 
at mwode af dfebam Iebtoauta seu ySnommos weeytsns sfbobust a 


Crs ae gm 


; 7 pa s+ >) 
10 »igae sesitq oft adnsadagex @ .sbuitiqms avey od atreassge® A aodw 


gies) ott .digasfovew 41d tinaeaiqes T bok, .aigia end mont ‘ide 
> 
bal lseriefavew navig o io ena wlgnta s19w golyog omits apag vatican . 
measfbotagq oft -ekaylens mexgubolzsq 02 19461 bas zasey tens ohnownid of 
Bbatiod aow 11 .esbudliqm: borhvpe sthodt Jeacteagn Incas aboiteq ie 
,ean7s2 ‘omld Inet to ytiolbolsaq to savanem sfdageny erev 8 od 09 
fenoboads nal) ese bos ,bInb ofa ot enokosutsull wobasx 03 seb yidiseog? 


gine land sonsiiay Yo revel mh 


ft) muitiosqe 5 26 , heselustng ot Yom 431488, g io holiettev out vag 
-lto20 Ic. 29fsnsipet? Io ywiatiasy © Yd ola Ssongi tev oil3 oF sok 2u8h 09 
jbellas of Yonoupay) qe yd muxtoege sd) of mordudf~atres edt “aovew 4 if 
-eblibe wid af Yondiiper? sits io soweq aia pera gD. soactmmas 


rwond al aalioe oma oft to mtiosgqa edd a30qmoOD oF hat 


a ae eee ane Ae 


A) IDEAL COS WAVE, PERIOD 277, PHASE ANGLE 0° 


WAVE LENGTH 


Y=cos 6 
WHERE OS @§ 277, 
AMPLITUDE 8 IN RADIANS 


BASIC INTERVAL 
27 


B) WAVE OF PERIOD 77, PHASE ANGLE 0° 
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C) WAVE WITH PHASE SHIFT @s5 


SINUSOIDAL MODELS FoR PERIODIC ANALYSES 


(AFTER RAYNER, 1971) 
FIGURE 2.05 
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as Variance spectrum, or power spectral, analysis. 


C. POWER SPECTRAL ANALYSIS 

In general, spectral analysis is used to examine a data series 
for evidence of predominance of certain scales of activity or influence 
over others in the underlying physical processes. The technique is 
useful in the determination of some order from within the usually noisy 
recesses of an often large and complicated array of data. Basic concepts 
in the method of spectral analysis and some fundamentals of the mathematics 
of the technique are presented here as background to the use of the 
technique in this research. More curious readers are referred to the 
works cited for details of the mathematical processes involved in the 
calculation of the spectrum. 

The notion of spectrum in spectral analysis implies the study 
of wavelength and period within time series. The spectrum estimates 
contributions made by both periodic and aperiodic functions within the 
stochastic series. It is initially assumed that the observed time 
series is a random sample of a continuous process which is composed of 
oscillations of all possible frequencies. A single series spectrum 
partitions the series variance into a number of bands or intervals of 
frequencies. The ultimate goal of spectral analysis is the determination 
of the nature of the generating processes such that predictions of future 
values, and occasionally so that simulations of historical records, can 
be achieved. A computation of the spectrum involves the calculation of 
the record mean, the autocorrelation function of the record, the Fourier 


cosine transform of each autocorrelation, and the estimation of spectral 
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density. 

he formation of the autocorrelation function is the basic 
operation in) spectral analysis.. In this operation, the series is correlated 
with avself using a variety of lags so that.a comparison of the fit_of 
the series to itself with its time origin shifted is obtained. Most 
analyses find a lag ratio (time shift of m units within a series of n 
points) of 10 to 20 per cent to be optimal for the analysis. That is, 
m/n < 20% seems the best general estimator of characteristics of the series. 
A. ten, per cent lag is used in all of the analyses conducted in this study. 

The measure of the autocorrelation is an ordinary linear 
coefficient varying from negative one (-1) to one (+1) through zero. 
The correlograms of Figure 2.02 again illustrate this method. 

Harmonic analysis is based on the classical theories of Fourier Series, 
which postulate that a series of trignometric generated curves, generally sine 
and cosine, may be used to represent any pericdic and transient functions by 
a series of trigonometric terms--generally as sine and cosine waves. Harmonic 
analysis may be defined as the determination of a finite sum of sine and 
cosine terms present in the record, and their relative contributions to the 
variation in the series. 

The "harmonic"! refers to the period observed when a specified portion 
of the data is considered. The fundamental or first harmonic has a period 
equal to the total length. The second harmonic will present one-half the 
period of the total, and the third will have one-third of the total period. 
Generally, the number of observations restricts the number of harmonics to 
n/2. Different harmonics are grouped or isolated according to frequency, 
and the importance of any particular frequency is then determined as its 


proportion of contribution to the total variance of the series. 
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The Fourier technique determines the harmonics by 


fitting sine and cosine curves, by 3 least squares 


method to the actual observed cycle. Generally, the first several 
harmonics describe .ne variation of the periodic function sufficiently 
well for purposes of interpretation of data patterns. 

The Fourier series may be represented as 


X fo] = 2 A [K]) cos (K 9 - @ [K]) where 6 is the phase angle, 


fe) 
or 360 t, for t representing any specific point selected from the total 


——_—s 


iE 
time T. The series expands as 


ajof{ cos O+ a tl (costh a a... Fa UR) -cocek 6 + o=6 
DlC\sineO bel] “sine 0. + b PR) “sin Ko ; 
thus represented as n/2 cosine terms and n/2 sine terms (after Rayner, 
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The variation contributed to the series by any specific harmonic may be 

2 

2 2 

expressed as __K_ , for Cy, = a ts be , and so C. has the amplitude of the 
Kth harmonic. The ratio of C, to the total variance explains the 


Proportional contribution of the specific harmonic. Since the harmonies 
are independent and not' correlated, no two harmonics can explain the same 
portion of the variance of X. Thus, contributions can be added so that 
af the firet harmonic contributes 20 per cent cf the total variance; 

the eecond, 00 per cent; and the third, 25 per cent, ‘then all other 
harmonics contribute only five per cent. 


The intensity of particular harmonics provides a measure of the 


42 


= at ee ner Milner 


*Tinary Ie « wt 
ibe Vo. nek 
picaany ain é 
stan Sendg on5 at Qiowedy cmt uw) ae: a 
terior ode wird Ge yuetew/ Janeg otha Rete Yee Nk: 
; 1 
seed 
0% fos. “Te * o.oo eae: tele 
OK he (ed © une 7 ae te Latea 
Toran eoyke) sioyed mnie Clin ban ems Satens ein ao boone | 
aa : fre wa 
: a 
af Ag yrange eaunsesi At 
1M 


Sas ie = : ae 
| (ak 96 
+! Z ~ } ts, yl pas ry ‘ont a J 


dy ie a Fronts hi! Ny enroya ANG es Sonat Fused a pains, 
old 26 Gbusi Jom of aes 40. be! Dery -. F Foy qo : 
oid GAlalons ssnaboiy [e263 ers Dn ra to bier otf a 
ehraeeted afi cst? «obtener 9) ¢ 2 54q- adi bender itn 


wikis ‘ls WlAlgxe oy aniiente os on frateleoyae Sen line 


(4oelawe iwigd ofr 82 dnvs-vog OD oo tint tain. sl Wena 
vale Lik ait , HS faq CE ,loeity ta bee qaden aay Oe" 


relative importance of each to the entire series variation. This may be 


shown in a periodogram, but the spectrum is a more versatile description 
of power since it decomposes aperiodic contributions, as well. 

Chance fluctuations may enter the series through the mode of 
recording, and will influence conclusions according to the difficulty 
bacountered in sampling techniques. Since these fluctuations could 
reduce the adequacy of the spectrum it may be desirable to minimize These 
effects. Various techniques or "filters" are available to amplify or 
reduce the apparent influence of certain frequencies within the data. 
Filtering is also used in detrending the original series. 

High pass filters allow the retention of only high frequencies 
in the data. Low pass filters similarly allow only the lower frequencies 
to remain in the series. Averaging is a form of low pass filtering as 
it removes high and some medium scale disturbance. Two common filters 
in spectral analysis are the differencing filter and the moving average 
filter. The differencing filter allows medium to high frequencies to 
pass unaltered, while lower frequencies are substantially attenuated. 

The differencing process may be represented as: 


oie (rye) (tL) 
where X eS is the value of the raw data at time t, and ree is the 
original value at the previous point in time, t-l. ee defines the 
residual data at position-t.~ Again, Cr Ey occurs in the residual 
analysis such that: 

Gig = 2 (em ty, *'O. Sb See) 
where £ is the specified frequency given in cycles per Qt (after Chang, 


1969). Moving average filters are similarly designed to pass high 


frequencies. 


’ i 8 a a9 [aS opie ain 
deen ad 
ad fa wf 6 rit Pent 3 int ad geass sad, ath) me ape. ty & Ti 1) ye ‘yD 
- - 7 
nel syisyoeob ols 168 aw ane & ai garsdaorgs aid sud ae 104 a a rod) 
vilow. a6 — bullies sale qs aseogmeasb 3P 2 x ¥Oq tc 


iy sbom of duis asitse oft t9ino yam amok sausoul?, 92 a. 


—_ 
volusiiiib 213. 02 asc SREP Ris aeraD gonou Tat ihtw le at sat (0397 


ane 


bivoyp enotjeusaeli saad eduhe -wavpindset gotiqnée eb ‘Somes 


Oo 


sts. 


“tod? ssimlatm oF sidagtesh ad yam JF muxtosge sty to east. ads oe 


~ e - 
o Viilans of oldatieve sts “e7ed Li", 16 asnpiniget Bice an 20339 
re re 
> 4 
' a Ke ee + ‘ F ¢ 
sRtab siz nidily aproneu petit Wregx95, 1° sarsusiint 7S eqD -. pry ’ 


seslive lantatrc eld aaa al wk bony oala ek 
- - Si : 
aoks at ote vireo Fodortnstbst Odo wells. 42 ed EE wis) fort vee 


eel si, ywol af} ylwo wolls ysralinta 954193 etanq wol ab fi a 
: Tas amt 
wired iti saeq wel To pfrol 5 af goignzvA sealjee sit ab nt aides, 
rie be ae 
2453117 nOmngp owl ssitglavsei) stlese murbem ames "bak be a 2ove 19% 2 
a 

~ 4 - * . Pp ; a ma 
TE19VH Bri von od 5re (otTis2 SB. 90HS7 tee e yt 932 ‘eure in Qa cf, 

t eytsisups tT? dat oF ampthem avelis «99529 gntone 1933 fh “at ioe 


-betevuesiis viletineiadue o1¢ asinasupor) +tswol olinw ‘neeediets 


i, sey - 
:e6 bedobastqed od Yoo eBeuddg pirtacera! ie | 
ia 


(r-7)% = aie “3 


ony wd X tog 7 onto 45 nteb wan oy Jo ou tay oe ae 


(I-29) 
sd4 ennideab st «[=3 ,9mr2 1) Inidg augive;q ada hal ow! 


_ igubtes+ Say at atwooo! (3) pb: thy, A e | on teg. ar 


: i 
: - - 3 a 


: hs 7 wet Gs Te & = 7 - 
7 Se 4) w= i > gas 
. ee Daan ie ae a Se ES 
- vs uae 19738) 104 dalaiio A ii movin 1 a0 Ups mga add ob 


el oe n << _ 
Snpglenpb views: > oh. 
7) Did ea pais di 7 eis 


_ 


Design of sampling interval and subsequent analysis must try to 
prevent "aliasing" in-the series. Figure 2.06 depicts the aliasing 
problem. The sampling interval equals or exceeds one half the period 
to be examined, and “hus the true period is never recognized. Only 


when the period exceeds twice the sampling interval can the event be 


measured, and the record length necessary to determine the event increases 


as the period of the eventapproaches a value of twice the sampling 
interval. The frequency (fn) which limits the events that can be deter- 
mined by any sampling frequency (fs) is commonly called the Nyquist or 
folding frequency, in which fn = fs/2. Sampling designs must consider 
this phenomenon carefully in order that all events considered plausiable 


are allowed a reasonable chance of detection, if possible. 


D. CROSS-SPECTRAL ANALYSIS 

Cross-spectral analysis extends the use of the technique to the 
simultaneous consideration of two series. "Correlations" of similar 
frequency components between the two series are tested by the "coherence" 
statistic, which is analogous to the squared correlation coefficient in 
correlation analysis. A plot of coherence against frequency can be 
useful in the identification of patterns of correlation, and therefore, 
of dependence between the two series. 

Cross-spectral analysis requires assumptions similar to those 
of its univariate analogue. Both series under consideration must be 
stationary, at least weakly, and random elements of each series should 
be Gaussian. Each series should be detrended before analysis. 


The cross-covariance function is analogous to autocovariance, and 
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——-—- ALIASED RECORD 
@ SAMPLING TIME 


TRUE RECORD 


THE PROBLEM OF ALIASING — OBSERVING PERIODIC 
EVENTS ViA SAMPLING INTERVALS OF LESS THAN 
ONE-HALF THE PERIOD (AFTER GUNNERSON) 

FIGURE 2.06 
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igs similarly basic to the cross-spectrum. For two time functions 
7 $3 ae hy Eat ies a eS 
Plt) and F2¢t) the cross-covariance is 
i \ T/2 
= T a) Eo : dt 
"XY (mn) =T/2 eyes Fern pee 


OF pit expected value notation, it is simply: 


Reveay =F LX ey HOM) > Ay] 


where m is the lag specification and iM, the series mean. 
The value of Ry is always a real number. Symmetry for the 
¥(m) 
cross~covariance occurs as: 


"EX CR): 


The eross-correlation function is obtained from the cross- 


covariance, and measures interactions between two processes. It may be 


represented as: 
R 
XY (mn) 


XY.) @===E 
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Its univariate analogue is the autocorrelation function. 

The Fourier transform of the cross-correlation of two series 
gives the cross-spectral density function for the series. The cross- 
spectrum is the complex product of the real amplitude spectrum (cospectrum) 
and the imaginary phase spectrum (quadrature spectrum). The cross- 


spectrum may be represented as 


to 


= 3 -2n 
*XY (us) 2 | Ryy (m) exp (-2Me mig.) 
oO 


with its components the cospectrum: 
©) 
eile CR + Roe 
XY (as) fs AE) YA(m) 


and the quadrature spectrum 
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The ratio ae.84 gives a measure of the phase relations between the 


Pyy 
two series, and arctan ( 7D) gives the phase angle Tear When coherence 
(£9 


estimates,arerhigh,orclose to.unity,. the, phase Pee) approximates well 

the relationship of the phase between the two series analyzed. However, 

low coherency is an indication of a poor relationship of phase and so 

a poor correlation at the specified frequency. (see Edge and Liu, 1970) 
Coherence is an analogue of correlation, but is given as a 


function of frequency as: 


C2 
XY 
Coh = (w) 
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It measures the correlation of the series within a small range of 
frequencies. Again, more detailed outlines of cross-spectral analysis 
may be found in Rayner, (1971), Matalas, (1966), and especially Jenkins 


and Watts (1966). 


E. APPLICATION 

Spectral analysis is a valuable technique in the development of 
the theory of river behavior. Spectra of both meandering and depth 
patterns allow some interesting comparisons, leading to theory, regarding 
the generating processes. Cross~-spectra of the two series provides 
valuable insight into a more encompassing set of processes. 

Spectral analysis is extremely useful in the condensation of 
masses of continuous data into meaning ful summaries of form, and is also 


important an that it can provide insight anto generating processes.. Its 
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use In Studies Of meandering is still in an exploratory stage, but it 
seems a promising technique, appropriate to the solution of many old 
problems. 

Chapter I described some of the work previously done in the 
examination of planform and depth patterns in rivers, and reference is 
now made to those studies which have used spectral techniques. 

Mackay (1963) first suggested the application of the spectrum 
to meanders, but Speight's (1965a, 1965b, 1965c) work was the first 
substantial applied contribution. Thakur and Scheidegger (1970), 
Chang and Toebes (1970), Yalin (1971), and Church (1972) all applied 
spectral techniques to meandering. Nordin (1971), Church (1970), Yalin 
(1971) and Chang and Toebes (1970) applied spectral methods to bed 
configuration downstream depth changes. The only work done so far 
utilizing cross-spectral analysis of meandering and depth is by 
Church (1972) and Chang (1969), although some of Yalin's theoretical 
predictions can easily be extended to cross-spectra. The findings of 
these works were summarized in chapter 1. Observations based on the 


data analyzed in this study, relevant to these works, will be presented 


in later chapters. 
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CHAPTER IIL 


COLLECTION OF DATA 


A. DATA NECESSARY FOR THE PURPOSES OF THIS STUDY 

The theory that the meander plan of rivers is best modelled by 
a sine-generated curve was formally proposed by Langbein and Leopold in 
1966. The theory postulates that a strong sinusoidal relationship 
exists between angular deviations of streamflow and distance travelled 
along the curved path of flow. It was also hypothesized that for rivers 
characterized by both meandering and straight reaches, the meandering 
reaches should exhibit lower variances in certain of their hydraulic 
factors (depth, velocity, local slope; also the planimetric geometry of 
the river). Langbein and Leopold (1966) assumed that the meandering 
river simultaneously maintains greater changes in bed contours than the 
straight river, the factors in combination finally rendering the meander- 
ing plan the more stable flow pattern. 

Meander planform has not commonly been examined over long 
reaches of streamflow. More usually, single meanders have been assumed 
to characterize the pattern (Leopold and Wolman, 1960; Langbein and 
Leopold, 1966), and thus single wavelengths have been chosen as sufficient 
description of plan morphology. Speight (1965) introduced both the 
study of long reaches of planform and the immediate observation that 
multiple wavelengths recurred to characterize the plan. Thakur and 
Scheidegger (1970) and Chang and Toebes (1970) have extended the analysis 


of planform as series of deviation angles, but data on changes of 
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directional angle of streamflow is scarce. 

Accurate information detailing the longitudinal profiles of 
river depths has been similarly scarce. Chang (1969) analyzed thalweg 
depths of the Wabash and White Rivers in Indiana using depth information 
on maps. Low-flow depths mapped at non-equal intervals were used to 
constrict a more continuous series of depth measurements through 
interpolation. Although it is possible to obtain continuous depth 
records of a river's longitudinal profile, no known depth analyses have 
yet exploited this potential. 

The paucity of adequate plan and depth information to corroborate 
or at least test such important theories as those mentioned earlier 
(Langbein and Leopold, 1966) emphasizes the need for further field 
study oriented towards a stochastic examination of the behavior of 
rivers. Particularly, analyses based on continuous depth information 
should provide useful comparisons involving both previous studies and 


existing theories. 


SB. CHOICE OF RIVERS FOR THIS STUDY 

Continuous depth records of longitudinal river profiles Cat 
high stage) and series of angular deviations of planform were obtained 
for three local accessible rivers (Beaver, Alta.; Kootenay, B.C.; South 
Nahanni, N.W.T.). The general locations of the rivers are shown on 
page 51 in figure 3.01 (a map). The rivers were chosen for their 
similarities to and for their differences from those of previous studies 
by other authors, allowing the best comparisons possible within practical 
limits. Co-operation of interested federal and provincial agencies, 


availability of adequate maps and air 
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photographs, existence of stage records and discharge statistics, and 
actual accessibility of rivers constituted the practical constraints, 
but these criteria fortunately interfered little with the choices made. 
The Beaver River exhibits regular meanders (within a spillway 
valley, relic of the Pleistocene glaciation) as it flows over a sand 
bed, thus providing potential for comparison with the work of Chang 
and Toebes (1970) on the Wabash and White Rivers. The South Nahanni 
River meanders, deeply incised, over a gravel bed, through very old 
canyons that have escaped recent glaciation(s). This river thus 
provides an interesting analogy for work on the Colorado River (Leopold, 
1969). The Kootenay River provides a useful contrast for these as it 
meanders irregularly, in some reaches entrenched, though not as deeply 
as the South Nahanni, and has a gravel bed. Detailed descriptions of 
these rivers will further emphasize their adequacies with respect to the 


purposes of this study. 


1) The Beaver River 


a) Geology: 

A meandering valley, cut by a much larger river active 
during ice melt, today contains and restricts the smaller (misfit) Beaver 
River (Neill, 1963). The old valley served as a spillway for large 
meltwater lakes. Its meander belt is about nine hundred feet (275 meters). 
The Beaver maintains a course of regular meandering within this valley, 
as shown by the map on page 53 (Figure 3.02) and again by the air photo 
in’Figure 3.06 on page 79. 


The Beaver River flows over a bed of clean medium to coarse sand 
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and between banks that ars generally silty but expose sand in places. 
Borings have shown the sand to continue for forty feet below the banks 
(Neill, 1963). The Beaver's average bankfull width of about 180 feet 

(or 55 meters) is maintained with little variation throughout the reach 
studied. Mean slope has been measured along the channel as approximately 
0.24 per 1000, and this seems to vary little with stage (Neill, 1963). 

An exceptional flood in 1962, during which the Beaver overtopped 
its banks by about ten feet, has been described as the Beaver's one- 
hundred year flood. The flood provided an excellent Opportunity for 
process study, and Neill (1963), who was primarily concerned is 
problem of pier scour, was able to observe some bedform evolution in the 
course of his study. A longitudinal sounding of a small portion of the 
river's midstream profile was obtained for near-bankfull stage. (At the 
time of the survey the Beaver was approximately one foot below bankfull 
stage and twelve feet below its flood peak). Two subsequent longitudinal 
soundings over a three-mile length of the same reach allowed comparative 
bedform examination. 

Neill (1963) observed, from the soundings, that shallows and 
deeps which had become very exaggerated during the various stages of flood 
had adjusted to form a less extreme profile by the time the Beaver 
returned to its normal low stage. Bedforms of three distinct types, 
classified by size, were recorded. Long waves (crest to crest approximately 
three hundred feet or 90 meters), dunes (about ten feet or three meters 
between crests) and ripples (dimensions in inches or centimeters) were 
observed, but these were regarded only as interesting observations and 
were not the main concern of Neill's study. 


Neill also recorded observations on some other flood effects. 
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sand were recorded on the insides of all sharp bends, 
while little bank erosion had been noticed. One bend had undergone 
considerable erosion on its right outer bank and simultaneous deposition 
on its left inner bank, resulting in a net conservation of river width 
(Neill, 1963). It is geomorphologically important to note both the 
occurrence of the flood and any possible resultant changes in planform 
(the only reference available notes that these were not substantial), 
because the maps on which planform data for this study are based were 
compiled from pre-flood air photographs, even though they were the most 
recent available. 

Large boulders are exposed along the Beaver's banks wherever the 
river attempts to cut laterally into the old valley fill of a past 
glaciation. The modern Beaver does not possess the competence necessary 
to move these rocks, which thus serve to protect or armour the river's 
banks. It is very likely that the boulders passively helped conserve 
the Beaver's planform during the 1962 flood. This would explain why 
almost no changes in planform were observed (Neill, 1963) throughout 
the various stages of flood. Such an armouring effect of these boulders 
would ensure the adequacy of the maps used in obtaining planimetric 


Gata. 


b) Vegetation: 

Low forest and dense bush cover the valley flats where cultivation 
is not widespread. A number of new roads into the valley have recently 
been constructed, however, indicating a trend towards expansion of valley 
cultivation. At the time of the river survey conducted for this study, 


most of the land adjacent the river supported willow (Salyx_ spp.), 
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aspen (Populus spp.), dogwood (Cornus stolonifera), and a variety of 


sedges and grasses. 


c) Climate 

Table 3.01 summarizes climatic data for the areas of the three 
rivers studied. Mean monthly precipitation and temperature figures 
for the Beaver indicate that spring runoff normally occurs during late 
April, continuing through May, while precipitation is heaviest through 
the summer months. Highest flow conditions are thus expected to 
characterize the stage of the Beaver through mid-to late May. The river 
will also flow at high stage (more sporadically) many times through the 
summer months (after storms). Neill et al (1970) have found that, normally, 
fourteen per cent of the Beaver's annual runoff occurs in April and 
another nineteen per cent occurs in May. Figures for runoff are similar 
for June and July, but in these months, runoff is likely more flashy. 


The survey of the Beaver was thus planned for mid-May. 


d) Flow Information 

The Beaver River flows northwest to join the Churchill River, 
eventually emptying into Hudson's Bay. It drains an area of 5460 square 
miles (14,100 taica: A manual guage, read daily, has been used since 
October 61955, toscolbect data on sthe flow ofethewBeaver River. Its 
location, near the Cold Lake Reserve, is marked on the map of Figure 
3.02, on page 58. 

The mean annual discharge of the Beaver River is about 936 cfs. 
(or 26.5 are}, Extremes have occurred on June 13, 1962 (the flood peak 


mentioned earlier) when a maximum instantaneous discharge of 21,800 cfs. 


Cor 617 ein was recorded, and on February 15, 1964, when a minimum daily 
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discharge of 13.0 cfs. (or 0.4 ma} was recorded. The mean annual 
floed flaws at 5,870 cfs. (or 166 ay Sediment data is not 
available for the Beaver River. 

During the time over which the survey was taken for this study 
(May 15, 1973), the mean discharge was approximately 962.0 cfs. (or 


3-1 


27.2 mS), slightly in excess of the annual mean flow. 


2) The South Nahanni River 


a) Geology 

Ford (1972) describes the South Nahanni Valley as the "most 
'fluvial' landscape in Canada". The region (Ford, 1972) escaped the 
Wisconsin glaciation so that during the past 300,000 years, fluvial 
processes have been the foremost processes acting to modify the land- 
scape. Isostatic rebound from the most recent glaciation, the “Parst 
Canyon Glaciation," has been responsible for the fluvial downcutting 
which extended the entrenchment of the South Nahanni River. Figure 
3.03 on page 61 presents a map of the river. 

The South Nahanni River meanders fairly regularly through four 
deep canyons within the southern MacKenzie Mountains. Small glaciers 
maintain the headwaters of the river. The fourth canyon terminates 
at its upstream end in Virginia Falls. The falls drop 300 feet (or 90 
meters) to the canyon's base. The recession of the falls, duringvtne 
past 300,000 years, is postulated to have been a significant factor 
in the development of the fourth canyon (Ford, 1972). Between canyons, 
the South Nahanni emerges as a braided river, but within canyon Gon Laine= 


ments it is very sinuous and meanders nicely. 
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immediately downstream from the fourth canyon, within oe First 
Dirarded SéctLron But at a pedrock constriction, the flow path turns more 
than 90°--so sharply that the flow diverges, forming a set of "figure- 
eight" rapids, with large standing waves present at the bend. However, 
throughout most of the remainder of the South Nahanni's course, the flow 
is fast but sufficiently deep to prevent the formation of major rapids. 
Figure 3.14 reproduces a photo of the figure-eight rapids. 

The depth of the canyon reaches varies between a few hundred 
feet and 3000 feet (or, between about one hundred and 900 meters). 
First Ganyon is the local base level for the other canyons, and it is 
the deepest of the four canyons (Ford, 1970). Maximum relief within 
First Canyon is 3000 feet (or 900 meters). 

The First Canyon is incised in a crystalline limestone formation 
(the Nahanni Formation) about 800 feet (or about 240 meters) thick, which 
in places overlain by Simpson Formation Shales. The limestone is the 
main cliff-forming unit in the canyon (Ford, 1970). There are no 
obvious glacial features within the canyon. First Canyon meanders 
fairly regularly, as shown in Figure 3.10, with a sinuosity of approxi- 
mately 1.5--the lower limit of "meandering" according to Leopold and 
Wolman (1960). Ford (1970) calculated the rate of canyon entrenchment, 
on the basis of dating of cave features and maximum relief, at no more 
than 2.66 feet (or, 0.8 meters) per thousand years, and he suggested 
that post-glacial uplift has likely accounted for at least 1600 feet 
(or 480 meters) of the canyon's relief. The gradient between Virginia 


Falls and First Canyon is very steep--about 2.25 per 1000. 
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b) Vegetation 
The South Nahanni River region supports typical boreal forest 
stands. Atop its cliffs and on any islands that persist for sufficient 
lengths of time, the soil supports mainly spruce (Picea spp.), larch 
(Laryx spp.), juniper (Juniperus spp.), wild rose (Rosa asicularis), and 


small plants such as bear-berry (Arctostaphylus uva-ursi), as well as a 


great variety of moss and lichens. 


c) Climate 

Table 3.01 lists monthly precipitation and temperature figures 
for the South Nahanni region. Maximum precipitation occurs between May 
and September, and temperature figures indicate that this adds to spring 
runoff in May and June. The river is normally frozen between October and 
early May. Climatic data suggest that the South Nahanni River maintains 
a high flow stage in early to mid-June, and so the survey of the South 


Nahanni was planned for mid-June. 


d) Flow Information 

The South Nahanni River flows south and east to join the Liard 
River which then follows a northeast path to join the north~- flowing 
Mackenzie River, eventually emptying into the Beaufort Sea. The South 
Nahanni River drains an area of 12,900 square miles (or 33,000 Km. ). 

A continuously recording guage, operating since October, 1959, 
collects data on the flow of the river near Clausen Creek. Its location 
is marked in Figure 3.03 on page 6l. 

The mean annual discharge of the South Nahanni River is about 
14,700 cfs. (or about 416 cu. m. per sec.) at the mouth of the First 


Canyon. (Surface Water Data, Water Survey of Ganada) This average 
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discharge is very similar to that of the Colorado River in the Grand 
Canyon. Mean monthly discharge for June, the month of highest flow 
stage, is about 62,000 cfs. (or, about 1755 cu. m. per sec.). Extremes 
have occurred on June 17 or 18, 1962, when a maximum instantaneous 
discharge of roughly 290,000 cfs. (or 8210 cu. m. per sec.) flowed in 
the river, and on March 25, 1961, when only 1690 cfs. (or, 45. ct. m. per 
sec.) was recorded as the minimum daily discharge. Accurate sediment 
data is not available for the South Nahanni River, but flow is very 
turbid throughout the summer months, and a boulder bedload is observed 
to move through the canyons (Ford, 1970). 

During the period of survey for this study (June 21 to 23, 1973) 
the mean discharge was approximately 24,798 cfs. (or 700 cu. m. per sec.), 


which is well above the mean annual figure. 


3. The Kootenay River 


a) Geology 

The Kootenay River follows an irregularly meandering path semi- 
incised between the Stanford Range and the Rocky Mountain Trench. The 
map given in Figure 3.04 covers the reach of the Kootenay studied for 
this thesis. 

Henderson (1954) published a detailed analysis of the geology of 
the Kootenay River region. For the convenience of the reader, the author's 
observations during the depth survey are supplemented with relevant 
information found in Henderson (1954). 

The upper Kootenay River is entrenched in a broad, deep valley 


that follows a major fault zone southeast to a point just downstream from 
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KOOTENAY RIVER 


FIGURE 3.04 


its confluence with the White River The valley in which the ‘upper 
Kootenay is incised maintains a fairly constant width of 2 miles at an 
elevation of 3250 feet. The river flows 50-300 feet below this main 
terrace level in its upper reach. Bedrock is visible in the banks. The 
upper Kootenay is, then, structurally controlled. Its general course 

is controlled by the occurrence of a narrow band of "soft calcareous 
chloritic phyllite" adjacent to the fault zone. (Henderson, 1954) 

Below its confluence with the White River the Kootenay turns 
to flow southwest, crossing another structural trend obliquely. In 
this reach, the Kootenay forms a sharp canyon at Gibraltar Rock. 

Further downstream, the Kootenay again changes its flow direction 
and proceeds south. The valley in this section becomes wider as the river 
passes over underlying gypsum. The new increased valley width is main- 
tained until the Kootneay reaches the flats separating the southern 
extent of the Stanford range from the Hughes Range to the east. About 
6 miles northeast of Canal Flats, the river tends slightly more towards 
the southwest and begins a braiding pattern which it continues much 
farther downstream. The river maintains its southwest flow direction 
adjacent an assumed oblique fault (Henderson, 1954) until it reaches 
the Rocky Mountain Trench. 

Large outcroppings of gypsum (sedimentary rock gypsum--Henderson, 
1954) occur along the Kootenay. Gypsum along the river reaches an 
exposed thickness of about 600 feet, and is interbedded with calcareous 
rocks of the Burnais formation (mainly limestone). Hoodoos occur along 
both banks of the river. 


There is evidence of Pleistocene glaciation of the valley. Thick 


deposits of stratified silts and gravels occur along the valley, overlain 
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by till in some places. However, the till is generally removed. 

Henderson (1954) suggests that the absence of u-shaped valleys, cirques, 
and hanging valleys in the Stanford Range seem to suggest that the area 
had been covered by a continental ice sheet that worked little to erode 


the rock. 


b) Vegetation 
Banks of the Kootenay are heavily forested, and extensive 
commercial lumbering is conducted in the region. Nearly pure stands of 


pine typify the regional forest. 


c) Climate 
Mean precipitation and temperature data for the Kootenay 

River area as,listed in table 3.01, indicatesthat spring runoff normally 
begins sometime late in April and continues through May. Precipitation 
is relatively high during May, and high stage is thus expected to occur 
in late May or early June. Summer storms should cause increased flow 
for short periods, but late spring melt high in the mountain headwaters 
of the Kootenay River is the main peak stage regulator. The survey of 


the Kootenay was planned for the end of May and earliest part of June. 


d) Flow Information 

The Kootenay River flows south to join the Columbia River, 
eventually emptying into the Pacific Ocean. It drains an area of 2080 
square miles. (or 5370 eee A manual guage (read daily) near Canal 
Flats has been used intermittently from 1939 to 1943, from 1944 to 1950, 


and from 1963 to the present. Its location is marked in Figure 3.04 


fy va 7 ye 


y2 Al 


eel ag m4 & AG ‘seen! 7 


show 03 +1298 hyttow rad avai mak Tetwintsong © 
—_ 7 : 
a 


XX 7 
ai “av 
oviurotke tae boteoto? ghivesd ots ysnsvoo! off To edna 
Yo abnétn s1ag vfs .cofysr add mt betovhass at gabvedmed fate 
sbi tani iq 
qerte sand act 10? 1560 aulistsgnes baa dotiesighoesg asat 
7 . ; > 
vlinwson PLones gnerqe Sadtestsstbat dO.¢ sidat at bevteti,da eeza, 


_ 
- 


notdmitgisert .geM teuerds acumtnoo bos fbaqgA at ste) ombsemee 
yoo 93 bedaeges buds. wt ogade dpit bow .yeM gabawb dgid yfoviss lor 
wisi peenoson) eeuso Nivota smote romwM@ .sowl ylass go yal sjel 
ar4ajavbeod alasdvom of9 nd dgid tiem gnisqe ete! tod ,abotisq trode 
79 ¢eveue sdT .voistuges sgaca Aenea alien od3 al sowid ysastood adZ 7 


nul Ins teeq teetitee bts ve 26. bno ari 202 bengalg sow yensdo0! & 


eolgservolal wolt | 

,1svi sidaulod af3 otol ot diuve ewok? reed yenatdood od? 

ORUS Yo ware mu aniedh 1T «96950 o) ibost odd Bye guierqes mi 
tenis sete Cylinh beer) sneug Fauitown A cm ore = asila 7 

Geet of MOLE avin? £h0L 09 C0! wes) «/ieoriertstat bow wiv “ 

| 09.6 wsvatt wl hedyamed noi isso) ast 1a oa ‘6 eee 


7 


Ts 


on page 66. 

The mean annual discharge of the Kootenay River is about 3140 
efs.\(or 90 cus m. per secs). Extremes have been recorded on May 24, 
1948, when the discharge was 29,700 cfs. for 840 cu. m. per sec.) and 
on April 1, 1941, when the discharge was a minimum of 178 Cis. or 5 cit, 
Mm. per sé@c.)./ » Themean: annualy flood: occurs in late May or early June 
carrying about 19,000 cfs. (or 540 cu. m. per sec.). Bed materials for 
the reach studied were observed to be gravel. 

During the time over which the survey was conducted for this 
study (May 31, 1973), the mean discharge was approximately 10,900 cfs. 


Govk300scashme perdisec.)eulThis isewell! aboveythe mean annual figure 


although it does not likely represent an annual flood. 


G. DATA GATHERED - @ SERIES 


1. Method of Measurement 

In order to compare planform geometry with theories that suggest 
that sine generated curves, sine waves, circular arcs, or other geometric 
lines best model the meandering habit of rivers, angular deviations from 
the path of flow must be recorded at constant intervals along the path of 
flow (interval denoted Ax). A chord is thus constructed between each 
point, Xe, and its adjacent point, Xeays located x units downstream for 
the series of n points. In this way, n-1 chords are constructed for each 
river. The angles formed between adjacent chords are measured and these 
are listed in order, forming the series of deviation angles. Thus each 


deviation angle, ¢@, is a measure of the angle between the chord joining a 


point Xs and its adjacent downstream point, Xipys and the chord joining 
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the same point x; and its adjacent point upstream, Kae Pasune 3.20. 
on page 95, explains the measurement of deviation angles in greater 
detail. Each Q.; then, depends on the location of the preceding point 
for its measurement, but is independent of any fixed origin or baseline 
for the series. 

Some interesting and apparently previously unrealized consequences 
of this manner of obtaining meander data will be discussed in Chapter 5, 
but it can be pointed out here that this is the method used by Speight 
(1965), Chang and Toebes (1970), Surkan and Van Kan (1969), Thakur and 
Scheidegger (1970), and Church (1972). These papers constitute all 
work published to date applying spectral analysis to the planimetric 
geometry of meandering rivers. 

This method of measurement for deviation angle of flow is an 
ideal parameter of meander planform for many reasons. The parameter 


¢- (angular deviation in radians) can be directly related to the radius 


of curvature R (derived as in Figure 3.20), for any segment of the river, 


as: 
ay PROBS 
R= 2 Sain (O-) 


This R can be used directly in the calculation of the centrifugal 


aeceleration,.)a,00f«the flow as: 


at 
aes R C 
where V is the flow velocity. Flow velocities were not recorded 


in this study, and so this concept has not been brought to use in this 
thesis. 
The use of flow path rather than specific origins or base 


lines (ie: down-valley Cartesian co-ordinates) in the calculation of 
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Q. eliminates the old problem of the possibility of multi-valued (non- 
functional) relationships. At the same time, it overcomes problems often 
encountered in the determination of a single flow direction or valley 
orientation to describe a river's path. A general valley orientation 
or baseline often becomes meaningless with respect to local flow 
direction, as the river's flow path can become complicated for tortuous 
or even very sinuous reaches. 

An arbitrary convention for this work designates left turns 
relative to the downstream direction along a river as negative curves, 
and right turns are considered positive curves. For convenience in 
later analyses, when removal of the mean is necessary for further 
calculations, right turn values (+4;) may be added to an arbitrary mean 
Value, wand leit tune values (-9,) may be substracted from a convenient 
mean, establishing a new non-negative series of D; valués. in this study, 
in order to facilitate computer coding, means chosen were 100 for the 
Beaver and South Nahanni Rivers, and 200 for the Kootenay River. Thus 
in illustration, for the Beaver, a value of 102 indicates a 2° turn to 
the right in a downstream direction, while a value of 95 indicates a 


O 7 
2 “turn to the left, as measured between adjacent chords. 


2. Thalweg Problems 


The river thalweg is the ideal line on which to base the measure- 
ment of both stream depths and planform angles. The thalweg traces the 
river's own main channel (and only channel at extremely low discharges) 
and so characterizes the satatnetae ability of rivers with respect to 


their plan and depth. However, information on the precise location of 


the thalweg is not available for many rivers. Few maps contain such 
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data, and few air photographs provide accurate information or.deta iling 
to expose the location of the thalweg. Fortunately, the midstream 
flow path adequately approximates the thalweg line, within the limits 
of accuracy allowed by most map scales, and so can be substituted for 
thalweg in the calculation of f;. That is, the difference between the 
locations of the thalweg and midstream paths introduces little error at 
the scale of measurement possible when common maps Or air photos are 
used. 

The problem of thalweg location is less difficult in the Pela. 
and so is more easily solved for the depth series. However, for large 
rivers, safe navigation of the thalweg during depth recording may often 
not be possible within the limits set either by the required accuracy of 


depth data or the navigator of the boat. Thalweg depth traces for this 


study were only gathered for the Beaver River. 


3. Maps and Air Photographs 


Adequately scaled maps and air photographs were available for the 
three rivers studied. All maps used were NTS maps, and all maps and 
air photographs were made available by the Federal Department of Energy, 
Mines, and Resources. 

The survey of the Beaver River, involving an 18.5 path-mile reach 
from the Alberta-Saskatchewan border to the Beaver Crossing Bridge, 
required only one map at the 1:50,000 scale. NTS map 73 L/8, compiled 
in 1968 from air photos taken in 1960, was found to be sufficiently 
accurate for the compilation of both the meander and the depth series for 
this study. Figure 3.02 gives the over-all view of the reach, including 


detailed plans and position fixes. 
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Maps of the South Nahanni were not available at the 1:50,000 
scale. The best maps available were the NTS quadrangles of Virginia 
Falls (95F) and Sibbeston Lake (95G) at the scale of 13250,000. These 
provide only crude location references. However, air photos of the 
entire reach were available at 1:1320 and these were adequate for loca- 
tional information. Photos, taken in 1956, are listed in table 3.03. 
The survey included 123.1 miles of the South Nahanni River. Figure 
3.03 on page 61 gives the overall view of the reach, including 
detailed plans and position fixes. 

Maps were adequate for the survey of the Kootenay River. 
General location was provided by NTS 823 at 1:250,000, while position 
fixes along the river were made using the 1:50,000 NTS series including 
82 J/1l2 W, 82 J/12E, 82 J/5E, 82 J/4E, and 82J/4W. The survey included 
46.1 miles of river from Canal Flats bridge to the “cable car crossing" 
upstream (approximately 50° 4S" O7"N, 115° 4S! O4"W). Figure 3.04 
gives the general view of the reach, while Figure 3.05 gives detailed 
plan and position fixes on the river. 

In all cases, for the planform series, x was chosen as one- 
tenth mile. The interval seemed to adequately cover river meanders 
while it simplified calculations. In retrospect, it probably would 
have been advisable to choose x as a function of river width, and this 
will be done in future research. The midstream line was chosen as the 
best estimate of river path, and deviation angles were thus found at 


intervals of 0.1 miles along each river's centerline. 
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D. DATA GATHERED: DEPTH SERIES 

Continuous records of each river's longitudinal profile were 
obtained for the purposes of this study. Continuous depth echo- 
sounding recorders were mounted in the boats, which then proceeded 
along the river on a path that maximized potential utility of output 
in relation to problems associated with each survey. 

A sixteen foot Larsen fiberglass boat, supplied and operated 
by the Research Council of Alberta,was used to survey the Beaver River. 
The boat was equipped with a side-mounted transducer which sounded 
depths with a beam angle of Le The recorder used in the survey was 
a DE 719 Raytheon, and a chart speed of two inches per minute was 
maintained throughout the survey. The original depth chart for the 18.5 
mile survey (series D) was 29.2 feet in length. 

Although a completely vertical sounding would have been ideal 
for the survey, it) was felt that the error incurred using the slightly 
foreward-angled transducer was insignificant for the purposes of this 
study. A thalweg depth trace was obtained for the Beaver River. The 
survey was conducted on June 22, 1973, approximately between 1500h and 
1900 h Mountain Daylight Time. 

Twenty-eight position fixes were recorded along the upstream 
path followed, and these were marked on the topographic maps and depth 
charts, but it was later found that two fixes were not accurate. Both 
fix problems were attributed to human error in the identification Or 
landmarks, and when omitted, did not seriously affect the record 
obtained. 

Figure 3.06 on page 79 shows, for a meandering section of the 


Beaver, the actual photographs used in this study, the depth soundings 
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Figure 3.06 Beaver River: section of typical regular meanders. Reach 
marked on air photo (top) extends upstream from mile 5.4 to mile 6.7. 
Chart (middle) shows longitudinal depth profile recorded as boat moved 
upstream. Cal comp plots (bottom) give digitized depth series adjusted 
to run's standard boat speed. 
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obtained, the position fixes used, and the velocity-corrected record 
for the same blocks. The correction procedure is outlined later. 
Figure 3.07 on page 81 illustrates the same qualities for a straight 
section of the rivei. 

Reconnaissance work on the Kootenay River indicated the necessity 
of a larger boat, and so a twenty-six foot Dorian jet boat, again 
supplied and operated by the Research Council of Alberta (constructed 
by Northwest Jetdrive in Edmonton), driven by a 455 h.p. Oldsmobile 
engine, was used for the Kootenay survey. The jet boat was equipped 
with both a Raytheon DE 719 continuous depth recorder attached to a hull 
mounted transducer (which again sounded depths with a beam angle of 
15) and a digital output second depth meter, which was useful in 
indicating the thalweg to the pilot. The Kootenay River survey was 
conducted on May 31, 1973, approximately between 1000 h and 1600 h 
Mountain Daylight Time (see Figure 3.19). 

The chart speed was again set at two inches per minute, and 
about forty feet of chart were obtained in the upstream run (Series CG). 
Thirty-seven position fixes were marked, and all were later used. 

Unfortunately, the Kootenay River was found to be too rough to 
permit a thalweg survey while maintaining safe operating conditions. 

A record of midstream depths was obtained instead (see Figure 3.18). 

Figure 3.08 on page 82 shows the air photograph of a curving 
section of the Kootenay. The original depth record for the blocks 
indicated and the boat-velocity corrected depth record are also shown. 
For contrast, Figure 3.09 on page 83 shows a fairly straight, narrow 
incised rock section of the Kootenay. Original and corrected depth 


records are included for the straight section, as indicated. 
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Figure 3.07 Beaver River: straight section. Reach marked on air 
photo (top) extends upstream from mile 11.7 to mile 12.4. Chart 
shows longitudinal depth profile recorded as boat moved upstream 
(series D). Cal comp plots (bottom) give digitized depth series 


adjusted to run's standard boat speed. 
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Figure 3,08 Kootenay River: entrenched section with some meandering. 
Reach marked on air photo (top) extends upstream from mile 16,2 to mi. 
18.4. Chart (middle) shows longitudinal depth profile as boat moved 
upstream (series C), Cal comp plots give digitized depth series for 
reach, as adjusted to standard boat speed for run, 
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Figure 3,09 Kootenay River: entrenched section with straight rock 
canyon, Reach marked on air photo (top) extends upstream from mile 
26.5 to 29,2, Chart (middle) shows longi tudinal depth profile, re- 
corded as boat moved upstream, Cal comp plots (bottom) give depths 
for reach, adjusted to run's standard boat speed, 
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Records were obtained for the South Nahanni River both as the 
boat proceeded upstream and downstream (series A and B, respectively). 
For the Nahanni surveys, a twenty-six foot jetboat was rented from the 
South Nahanni River Expeditions Company, based in Fort Simpson, N.W.T. 

A tube type Raytheon depth recorder was used, supplied by Water Survey 
of Canada, with a side-mounted transducer, again having a beam width of 
DST eka upstream survey of the 127 mile reach began June 21, 1973, at 
approximately 1600 h, continuing until about 2200 h the same day. 

On the second day, the survey began about 1000 h and continued in 
heavy rain all day. A downstream survey was conducted on June 23, 
beginning about 1000 h and continuing until about 1900 h. Chart speed 
for both upstream and downstream surveys was set at 36/60 inches per 
minute. The downstream chart produced was only 4 feet long, and over 13 
feet of chart were recorded during the upstream run (see Figures 3.16 
and 3.17). 

Position fixes were taken for both the upstream and downstream 
runs, in identical places where possible, and a total of eighty-four 
were recorded during the longer upstream survey. A number of the fixes 
recorded in the braided sections were later found to be unusable. Island 
changes during the time since the original production of the air photo- 
graphs confused the two navigators. Landmarks were not identifiable from 
photos, and there was no immediate solution to the problem. Consequently, 
individual data "blocks'' are very long for the braided section. 

A thalweg survey was considered unfeasible for the South Nahanni 
River. The boat operator's original opinion was that the thalweg would 
be too swift to allow safe collection of data. In retrospect, however, 


it seems that a thalweg trace could probably have been obtained for the 
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river, and its is unfortunate that this was not done. Instead, the 
river's “anti-thalweg" was followed in order that quieter water might 
be navigated. Figure 3.15 exemplifies the rough-water problem. 

Figure 3.10 on page 86 shows, for the meandering section of the 
First Canyon, the blocks used, depth chart obtained, and corrected 
depths. Figure 3.11 on page 87 provides similar information for the 
Second-Third Canyon areas. For the braided sections, Figure 3.12, on 
page 88 shows these features. (In this figure, block #10 appears 
elongated on the chart, as a slower chart speed was attempted for this 
block only.) "The Gate", a spectacular narrow canyon in Third Canyon 
is shown in Figure 3.13, on page 89. Figure 3.13 again gives both the 
original and the velocity-corrected depth record for this section of 
the downstream run. It was not possible to obtain a record for this 
reach in the upstream run as velocity problems overcame the depth 
recorder. 

The Digicon digitizer, built by the University of Alberta 
Technical Services Department, was used to digitize and file the depth 
chart data. In the digitizing process, the chart is placed on a 
magnetic digitizing table. A programmer then traces over the actual 
record with a moving arm that magnetically records its position to the 
nearest specified interval. Charts in this study were digitized at a 
Ax of one-fiftieth of an inch. It was felt that no usable information 
was present in smaller wavelengths than this sampling interval provided. 
Individual blocks (the chart lengths between subsequent fixes) were 
digitized separately and stored in separate files on a single magnetic 
tape. The digitized record was transferred to an IBM compatible 9 track 


magnetic tape for permanent storage. 
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Figure 3.10 "First Canyon" of the South Nahanni River, Reach marked 
on air photo (top) extends upstream from mile 82.2 to mile 76.9 below 
Virginia Falls. Chart (middle) shows longitudinal depth profile for 
reach, recorded as boat moved upstream (series A). Cal comp plots at 
bottom give digitized depth series adjusted to standard boat speed. 
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Figure 3,11 Second and Third Canyons on the South Nahanni River, Reach 
marked on air photo (top) extends upstream from mile 33.3 to 27.4 
below Virginia Falls, Chart (middle) shows longitudinal depth profile 
recorded as boat advanced upstream (series A). Cal comp plots (bottom) 
give digitized depth series adjusted to standard boat speed for series, 
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Figure 3,12 South Nahanni River braided section, near Liard confluence, 
Reach marked on air photo (top) extends upstream from mi, 114,2 to 108,8 
below Virginia Falls, Chart (middle) shows longitudinal depth profile 
recorded as boat advanced upstream (series A). Cal comp plots (bottom) 
give digitized depth series adjusted to standard boat speed for run A, 
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Figure 3.13 "The Gate" in Third Canyon on South Nahanni River. Reach 
marked on air photo (top) extends downstream from mile 33,2 to 40.6 
below Virginia Falls, Chart (left) shows longitudinal depth profile 
recorded as boat moved downstream (series B). Cal comp plots (right) 
give digitized depth series adjusted to standard for river. 


E.” DATA’ GATHERED: CORRECTED DEPTH RECORD 

It was™necessary, in the initial stages of analysis, to keep 
individual blocks separate, as it was known that the boat speed must 
vary along the path travelled (both for logistic and safety reasons). 

As the boat proceeded upstream, shallow areas would force the boat to 
proceed with a forward velocity slow relative to its mean velocity. 
Deep pools in the path would cause the boat to move relatively quickly. 
The depths were recorded at a constant time interval by the depth 
recorder, and thus at higher boat speeds, the same length of river 
would produce less chart than would be obtained at slower speeds. It 
was necessary to correct for this error during the data analysis. 
Since the ratio of the river distance-travelled to the corresponding 
chart distance obtained could be found for each block (between two 
known fixes) the mean velocity travelled along any block's river distance 
distance travelled. 
could be easily computed as time elapsed . Elapsed time was found 
in terms of inches of chart comprising any block compared to the 
corresponding chart speed. 

Once the mean boat speeds were calculated for each block, the 
block corresponding to the fastest mean velocity in any survey could be 
determined. All other blocks in the survey then contained relatively too 
many points, compared to the "standard" block, and points could then be 
systematically removed from all other blocks to correct for velocity 
fluctuations. 

The number of points removed was proportional to the ratio 
between points per unit distance along the river in the standard (fastest) 
Dicck and tthat of the ae compared. It was felt that a systematic 


removal of points would not significantly modify the areas of interest 
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of the spectrum of the block involved, considering that the x of the 
chart record digitizing was one-fiftieth of an inch. Point removal 
conducted systematically only corrects the record, so that blocks can 
finally be joined ac a continuous series. Appendix 1 on page 185 
lists the original number of data points in each block and the corres- 
ponding percentage of these removed by this process. The removal 
process is outlined in further detail in the appendix. Figures 3.06 - 
3.13 allow visual comparison of many blocks before and after point 
removal. The figures show well the success of the method outlined, 
and the lack of “losst' of relevant data. An alternative method of 
correcting for boat speed fluctuation would be to "stretch" some (or 
even all but the slowest block) blocksusing some means of interpolation. 
Either (or both) method would probably give satisfactory results, but 
point removal is simpler and no significant information is lost if x 
is sufficiently small. 

Once the corrected continuous data series was obtained, tape files 
were corrected, and the final correct series were stored. This data may 
now be obtained upon request, from Ms. D. Dodd, Programme Analyst, 
Geography Department, University of Alberta. A blank magnetic tape should 
be supplied with each request. Data will then be duplicated onto the tape, 
and the tape will be returned as directed. Alternatively, the uncorrected 
data may also be obtained. The original chart records will be given to 
and stored by government agencies (Beaver: Research Council of Alberta: 


Edmonton; Kootenay and South Nahanni Rivers: Water Survey of Canada.) 
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Figure 3.14 


"Ficure-eight" Rapids on the South Nahanni R. 


Figure 3,115 


Rapids on the South Nahanni R. -- series A, 5 
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Figure 3.16 Transducer mount for 
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CHAPTER IV - ANALYSIS OF DATA 


STATISTICAL PROPERTIES OF THE SERIES 


AQMOSTATIONARITY .OF> THE SERIES 

The Pores of stationarity, at least in its weak sense, 
is necessary to analysis by variance spectrum statisticse This assumption, 
discussed earlier in Chapter 2, essentially provides that all properties 
of the series obtained remain constant throughout the ensemble. The 
sample thus shows independence from absolute time. A series that 
demonstrates trends in its means or variances is said to be non~stationary 
or evolutive. 

Stationarity in the weak sense implies that the series may be 
characterized by its mean and variance, which do not show significant 
change over the length of the record. Assumptions of higher-order 
stationarity, in which higher order moments are also constant throughout 
the sample, are rarely found in nature. 

Trends constituting non-stationarity of a time series may occur 
in the means or variances, or both. Simple filters quickly remove 
trends in the mean. However, trends in the variances along a series 
require more difficult manipulation. Deterministic elements must be 
calculated for the series and removed. 

Downstream increase of stream discharge and river dimensions, 
especially present if the series contains major confluences, diminish 


the probability of stationarity for the series. In this study, both 
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the South Nahanni and Kootenay River reaches included confluences of 
unknown significance with respect to this property. The Beaver River 
reach contained no substantial confluence, although the river exhibited 
a very slight tendency to increase in width and depth downstream. 

Each depth chart and planform series obtained for the study 
constituted a realization (sample function) of the stochastic process. 
Strict and even weak stationarity are statistical properties used to 
describe ensembles of time series. The only stationarity test appropriate 
for simple realizations of the ensemble is that of self-stationarity, 
which analyzes changes in properties along the sample series. A series 
is thus self-stationary if subseries derived from the main sample by 
simple segmentation demonstrate no significant changes in their moments 
when compared as units. Weak self-stationarity implies constancy, 
within random sampling fluctuations, of the first two moments. Yevjevich 
(1972) suggested that subseries with parameters confined within 95 per 
cent tolerance limits about the corresponding values of parameters for 
the entire realization, may be inferred to be self-stationary. 

Table 4.01 presents chi-squared (x7) values obtained for tests 
of self-stationarity of depth variances and of both relative (~ 9) and 
absolute (}9\) curvature series. All series (A, B, C, and D) were 
tested and all were found to be weakly self-stationary. 

Tests for self-stationarity of depths aataiied initially the 
segmentation of each digitized corrected depth series into 100 blocks 
of equal sizes. Means and standard deviations were then calculated for 
each block, forming four new series (A, B, C, and D) consisting of one 
hundred block-mean depths each, and four new series (AVEC. and UD) 


consisting of one hundred block-standard deviations each. Frequency 
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histograms were constructed for each of these series. These are shown 
asm teures: 4.01, 4.02, 4.035, 4.04, 4.05, 4.06, 4.07, and 4.08... The 
frequencies found were used as expected values in the x tests performed 
lager 

Each series was next considered as composite of two unit samples 
of equal lengths (50 blocks each) and these halves were then examined 
for similarity to each other and to their parent series. Histograms 
of each half-sample are included in figures 4.01 through 4.08. Frequencies 
found were considered observed frequencies for the xe? test. Expected 
frequencies were adjusted to correspond to the half-size of the samples 


2 , 
tested, andX values were found in the usual manner, as 


where e, is the frequency expected, 0; 


; is the frequency observed for each 


i'th interval, and n specifies the number of frequency classes constructed 
for the test. 

Statistical tables (Walpole, 1968) furnished critical test values 
at 95% tolerance and Y= n-l degrees of freedom. Degrees of freedom 
were calculated as n-1 since only the total number of observations was 
used in the calculation of expected values. The test thus ran: 

Hos fy = 45 = se or 81 = 8607 165 

Hy: 1 ao Fhe or OF 7 Oy 19 5 
where m4 and 6, are the mean and standard deviation of the first half of 
the parent series s, and M, and 6. are the parameters of the 
parent series. A value ote greater than the value tolerated for 95% 
of occurrences at n-l degrees of freedom would cause rejection of Ho 


as beyond the limits of chance. 
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All first and second order moments supported the assumption 
of self-stationarity of the depth series obtained, as seen in Table 4.01. 
All smaller reaches also exhibited weak self-stationarity at 957. 

Meander data were also tested re statistics. Zach 
observation was recorded atOx = 0.1 miles, and for the x7 test, it 
was sea He that each angular derivation was thus a mean for 0.1 miles 
of river path. Circular paths were assumed on the micro-scale used 
in thevcalculation of radius of curvature, and so variances of angular 
deviation fluctuations, with the 0.1 mile path, were assumed constant. 
Both relative and absolute curvature series were tested in their 
"means'', then, and all were found to be well within tolerance limits 
of 95%. Meander planform series were thus also found to be weakly 
self-stationary, as indicated in Table 4.01. Figures 4.09 through 
4.16 depict histograms of the meander planform data. 

Results obtained from spectral and cross-spectral analyses of 
of the data may thus be treated with greater confidence under the 
assumption of self-stationarity. The data thus allow valid use of the 


statistical technique chosen. 


B. INFORMATION CONTAINED IN THE RADIUS OF CURVATURE 

Radius of curvature was calculated, as shown in Figure 3.14, for 
all deviation angles measured for the planform analysis of the rivers 
studied. Calculation of this statistic was based on the assumption that 
for any three adjacent points in the planform series, a’circle fit to 
the points reasonably approximated the flow of the river. Figure 4.17 
presents graphically the radius of curvature statistics for the study. 


Radius of curvature was chosen in order to examine preferred 
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angular deviations, as given in table 4.02, in non-dimensional terms. 
Radius of curvature was divided by river width for each reach to provide 
a dimensionless value for study. It was found that for reaches with 
meanders or curves, the parameter R/W clustered in a small range of 
values. River size or scale factors thus seemed eliminated when the 
statistic was used. 

The parameter R/W showed an interesting tendency to approach the 
value 271as bend tightness increased (ie. as the sinusoidal form was 
approached). The significance of this observation is difficult to 
ascertain in light of the "approximate" nature of the values used. 

Figure 4.17 presents observations on the statistic R/W. It 
shows no distinct relationships with width or depth of the rivers, 
but size of sample (small) and accuracy in deviation of R/W make 
inferences difficult. The figure e included in order only to present 
the data as it was found. 

The values presented can only be regarded as rough approximations 
of their true values. Widths were obtained from maps and air photographs 
rather than in the field, and R was calculated for angular deviations 
accurate to one degree. Inaccuracy involved thus in the computation of 
R prevented its use in more sophisticated analysis, but it seems a good 
parameter, as it directly relates to flow acceleration in bends (centri- 
fugal acceleration, a, where a = w/e, and v is the velocity of flow.) 


Further work with the parameter is planned, but remains beyond the scope 
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G. INFORMATION CONTAINED TIN SIMPLE FIRST AND SECOND MOMENT STATISTICS 
Simple first and second moment statistics can be used to describe 
a limited range of form-process phenomena related to both the bed 
configuration and the planform geometry of rivers. The mean and 
standard deviation may be useful in distinguishing reaches within a 
longer series describing the river, where certain reaches are comprised 
of a number oy blocks, the characteristics of which are similar within 
each reach, but significantly different between reaches. ‘The information 
contained in the means and standard deviations of the depth and angular 
change series provides a general overview of changes characteristic of 


the river's behavior along the path of flow. 


Ibs Depth Series 


a) General: 

Means and standard deviations were computed for individual 
blocks in each "corrected" depth series (chapter 3 explains the correction 
performed on the raw data). Blocks were then arranged according to the 
different geographic reaches contained within the whole series, and 
mean depths and average deviations were calculated for each reach. 
Table 4.04 presents these statistics (first and second moments) for 
all the rivers. 

A strong pattern is evident in this data. Straight sections 
of channel are characteristically shallow relative to the river's mean 
depth, and when not entrenched, they exhibit substantially smaller depth 
variation along the flow path than other reaches. Sections with bends 
exhibit greater mean depths and depth variations, and regularly meandering 


sections show both the greatest mean depths and the greatest variations 
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Table 4.04:A FIRST AND SECOND MOMENTS OF DEPTHS FOR REACHES STUDIED 
FOR SOUTH NAHANNI RIVER ~- includes both series A - upstream, 
and series B - downstream 


Depth Statistics 


3 in feet 
Reach x Sed. Comments 
1. Series A - entire LO 782 Sree 
main braids pels) 6.98 A: greatest variation 
(mi./81-115) 
First Canyon L2.09 4.03 A: most deep 
(mi. 68-81) 
intermediate braids 9.12 geal 
(mi. 59-68) 
Second Canyon O29 ae 
(mi. 49-59) 
Third Canyon LOoo 4,20 
(mi. 27.5-49) 
flats 8.83 3.06 
Oma © 7-27.5) 
Fourth Canyon Ds 02 be A: most shallow 
(mi. 4.5-7) A: least variation 
2. Series B - entize™ 11.74 yar ial 
main braids 1 ae | eihe) Bs greatest variation 
First Canyon 12.49 4.82 
intermediate braids 10.08 ae bd B: most shallow 
Second Canyon 12.94 3.42 
Third Canyon 13.01 5.04 B: most deep 
flats LO. iS Zend B; least variation 


al 


“note, series B excludes Fourth Canyon and half the flats section, accounting 
for discrepancy with series A. 
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3. 


FOR BOTH SERIES C - KOOTENAY RIVER AND SERIES D - BEAVER 


RIVER 


Reach 


Series G - entire 


braids 
Cmieitl= 525) 


semi-entrenched 
(mi. Seo et) 


straights section 
ie 20 15-2505) 


irregular meanders 
(mi. 25.5-26.5) 


rock canyon 
(mi. 26.5-29) 


shallow- incised 
irregular meanders 
(mi. 29-47.1) 


Series D - entire 


meanders - border 


(mi. 1-9) 


gentle curves 
(mi. 9-9.7) 


straight 
(mi. 9.7-13) 


meanders 
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Depth Statistics 
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O20 


7.50 


D.5 


Soke 


Tat 


Dit Vill 


4.27 


aol 


in £Eeet 
Side 


2.77 


1.84 


1.80 


20D 


1.63 


0.332 


1.386 


Comments 


most shallow 


beast variuation 


most deep 
most variation 


most deep 


greatest variation 


most shallow 
least variation 
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in depth. In general, for non-incised reaches, depth variation tends 

to increase as tightness of bend increases (ie. as radius of curvature 
decreases), implying a relationship of greater variation in hed 
configuration for reaches of greater variation in angular deviation. 
Table 4.05 lists the rank order of depth variations and the corresponding 
rank order of angular deviation variations. Braided sections generally 
show greater depth variations than either the straight or irregularly 
meandering reaches, though their mean depths are less-than-average in 
_the context of the entire river. Entrenched sections, where flow is 
confined between rock banks, as in the cases of First, second, and Third 
Canyons of the South Nahanni, and the upper reaches of the Kootenay, 
also show both large variations in depth along channel and deep channels 


relative to other reaches in the same river. 


b) Beaver River: 

The portion of the Beaver River studied for this paper consists 
of two separate reaches of regular meanders, misfit within a larger 
valley, joined by more gentle curves of channel which continue to become 
a rather straight path of flow. Figures 3.06 and 3.07 include air 
photographs, actual depth chart, and depths corrected for boat velocity 
variations from these reaches. The meandering reaches exhibit greater 
depth variations than either the straight section or the sections of 
gentle curves, and the mean depths are correspondingly greater for the 
meandering reaches. The reaches of gentle curves are both deeper and 
more variable in their depths than the straighter reach. Downstream 
meanders tend to be generally deeper than their upstream counterparts. 


Variations in the depth profile of the meandering reaches tend 
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Table 4.05: RANK ORDERING OF VARIATIONS IN AND DEPTH FOR ALL REACHES 


(greater variance assigned lower magnitude rank) 


A. Series A and B, South Nahanni 


Reach Ald| AZ® A depths BJof Bt 8B depths 
1. Main braids 1 1 1 1 i t 
24 First Ganyon 3 3 2 iS 3! 3 
3. intermediate braids 2 2 5) ye 2 4 
4. Second Canyon 4 4 6 4 4 5 
a3.s third Canyon 4 4 4 4 4 2 
6. flats = - 5 - ~siateatdegs 
7. Fourth Canyon 5 5 7. 5 J es 


Be Senewacy [¢ 


Reach c || Ci . G depths 
1. braids 3 2 4 
2. below White River i 1 | 


3. above White River: 


- straight section 6 
- irregular meanders 2 
- rock canyon 2 3 1 
- irregular meanders 5 
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to be explained more by a regularly spaced pattern of well-separated 
pools and shallows, as may be seen on the depth chart in Figure 3.06, 
than by smaller scale, higher frequency bed undulations. Pools occur 
regularly at the sharpest portions of bends along the meandering sections, 
while shallows are more pronounced at points of inflection. In contrast, 
the depth series of the straighter section consists of higher frequency, 
lower amplitude oscillations. The Beaver is a sand bed river, and ge 
variance in the straight section is explained, at least in part, by the 
occurrence of dunes. The air photo in Figure 3.07 shows that dunes 
characterize the bed in straight portions of channel. Figure 4.18 allows 
quick comparison of the actual depth series for the contrasting sections, 


while listing corresponding depth statistics. 


c) Kootenay River: 

The portion of the Kootenay River studied for this paper begins 
in braids, but upstream the river incised into irregularly meandering 
sections broken by straight sections. The two main straight sections 
provide a good comparison of bed configuration, as one section is an 
incised rock canyon while the other is confined by a high rock wall on 
its left bank, but flows past a lower terraced right bank. The rock 
canyon was both the deepest and most depth~variable reach, while the 
semi-confined straight reach was among the most shallow and least 
depth variable for the river. The irregularly meandering sections were 
generally deeper than the river's average flow depth, and they showed 
more pronounced pool-and-riffle sequence in the river series. The braids 
form the shallowest section of river. 


The highest frequency of depth fluctuation (in long profile) 


enoTI2da grtashapem sia, sate ond ‘a: siesta 
wjuastvon af .noigoelhet Do aieloq ta beaneono7g stom Sta 
\yonwupaat saitgiit 2b" Selina, dakqoee titgtexde ‘ede So! eababa! 
aéy boe ,vovbu bad Apa w at xévael sdT senotsed It oe obu: it iqn savel 

eds vd ,iteq fi Jasol Ja ,baminiqus 23 nots 208 sdiezse of pt sonetaay ; 
anid tad% owors VO.E swugit mb orodg ate oft .ennb oe sanetauosa, — 
ewolle Bi.d s10g]% Joniads Ye andtrteq Sdglavse at bed odd osisasouxeda 
(ShOl6298 gta etIPS of) Ww setros aged caaiiiass aij to noattsqmos x>tup 


.aatverieie diqab: gribinog nea gniint! oitiw | 
» git ee 
‘i 9veN esnor00t Ca, 


eniiged is¢aq ghd? 103 beibyyea sovid yonesond adj to nold10g edt Saul 


_fetrhassm yiveluneys! otnt boatoat tovia sé3 meodzequ -Jud _eblead me 

evaltooe jdyteria nism owt oT .anoltoes: tdgtaxie a, astord enol t208, ' 

of ai MOltoog ano te ,polanwagliner bed To noarasqmos boog s abivory 

we Iiaw toon 4gif a wt banl@noo al rerdto ada oltdw nuyaks Agvet bostont 

jogos oft Jind Ivighs besexxe2 iswol « teeq ewold aud eined Jot 23% ; 

ods ofiiw . fons aidaispy-daqeb jaom bas jesqoeeb of taod caw HOLGER: 

jeee! bon wollpde Pear of2 grome esw donst idgtnuda bamiInos-kmae ' 

aay wrtoiiona irhieiicem qiaelogerst sdT peck any x08 mua jb 

bewody vats one , Adega wed Syeteve a! set aif ad pies nog, 
ebtaa! eff .évicve seve ed al soneupse o»J1T37La-bns-I 

-tevit Jo nedtooe 4 


cine anol. ak) aokinesou! digeb te qanauperd 38 


= 


i 
1. 
n ry in| 
1.) 


2 
= 


Sty Sang 


LHSIVYULS DSNIYVdNOD ‘SSN F=XV LY ONIIdNVS HIM ‘SHINES HLdad WNLOV 


SATW FF =XV “HLVd SNOW AONVILSIG YSsAlddNn 
8E 9E ve ce O€ 8¢ IE v2 ce 02 8 Sp vi Gil Ot 8 2 v @ 


eke 


Lassa SHidsG 


vs 


SHacinvyawn WVadisd(> ——— 
HOVse LHSVels 


8°8 


A@ BEVCHES 4 LHEIS BED COVLICNHYLIONS” BEYAEL BIAEU 


REMARL DIELVMCE WrOve Lyin’ Ux=‘t wirke 
/ 
SJ 


4d OM NOS 1% 


—> 7c 
fm, ei 
cil i Ral 
eae e+ 


ibe: 
i 


te 


s) 


ge i) ‘law “ 


a i 404 : : 
r 


—~ j iY Oem rT 


i @€ neh - 


35 


i 


a ; : . 
( ; 
7 . 7 1 ie 
a ee 
< +e Galas 
hg we . ad 
x ye ere 
“fs a : 
m aad } ais 


j ’ i” gives “ = 


a. boned 


occurs in the rock canyon, shown in Figure 3.09, where fluctuations 

are also of great amplitude. Flow is most turbulent in this section, 
where four-foot standing waves made upstream progress difficult even 

by jet-boat. Not far downstream, the less entrenched straight reach 
depth profile exhibits a lower frequency, lower amplitude series of depth 
fluctuations. Flow was less turbulent, and the river was wider here 

than where it became more entrenched upstream. 

Upstream from the rock canyon, the river remains entrenched, 
though not as deeply, and it begins a pattern of irregular meandering. 
Depths in this section, according to the profile obtained, vary with 
high amplitude and intermediate frequency. The deepest pool (44 feet 
or about 14 metres) was recorded in this section, while pools here 
averaged a depth of fifteen to twenty feet. The width was confined, 
though not canyoned, and flow was again rather rough, characterized by 
standing waves. 

The braided section begins where the Kootenay emerges from 
entrenchment and flows through a broad valley. Depths in this braided 
section are generally less than the river's mean depth, and bed undulations 


Occur irregularly in the bed. 


d) The South Nahanni River: 

The South Nahanni River, below Virginia Falls, consists of four 
deeply entrenched reaches joined by braided channels, finally opening 
into intense braiding that continues to the Liard confluence. Three of 
the canyons meander irregularly, while one is comprised of rather regular 
meanders (comparison provided in Figures 3.10 and 3.11). First Canyon, 


the most regularly meandering and most deeply entrenched section, exhibits 
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much greater-than-average depth and depth variation in the river 

series. The less regularly meandering Second and Third Canyons also 
possess very deep channels and highly depth-variable bed configurations. 
Third Canyon, in which The Gate is located (shown in Figure 3.13), shows 
both the greatest mean depth and greatest depth variation for the river. 
The braided sections show slightly less depth and depth variance, but 
are deeper and more variable in their depths than the irregularly 
meandering Fourth Canyon just below Virginia Falls. Fourth Canyon 

is the most shallow and least depth variable section among the reaches 
studied. It is the farthest upstream reach, and thus directs the 
smallest volume of water. 

First, Second, and Third Canyons vary in their depth series 
with large amplitude oscillations related strongly to the occurrence of 
bends. Sharper bends show deeper pools and again, the "shallows" 

(here, about ten to fifteen feet in depth) occur regularly at points 
of inflection. Higher frequency, lower amplitude fluctuations appear 
in the series, but are noise.and their peaks are not particularly 
related to position in the bends. Variance of the series is mainly 
controlled by a pool-and-riffle sequence. 

Braided reaches between canyons exhibit higher frequency, lower 
amplitude depth fluctuations, and show slightly higher depth variances, 
but smaller values for mean depth. ‘The braids below the canyons 
approximate the rivers mean depth. This reach is shallow relative to 
the canyons, and marginally less variable in depth-Undulations of 
intermediate frequency and amplitude show only little relationship to 


planform within the noisy data array. 
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e) Theoretical Implications of Depth Observations, Using Mean and 
Variance 

Observations based on the first two moments of the depth series 
data gathered for this study support the theories of Langbein and 
Leopold (1966) and Leopold and other (1960) with respect to differences 
in bed configurations of curved and straight reaches. Data discussed 
eee for rivers of small, intermediate, and large sizes, 
and a variety of channel patterns and bed materials were considered, 
though briefly. 

A pattern of increasing mean depth and depth variance with 
increasing regularity of meandering and increasing intensity of entrench- 
ment is definitely established within the depth data. The increasing 
depth and depth variance patterns shown are consistent with the theories 
of equalization of energy dissipation through unit lengths of river, 
which propose that curves introduce additional organization into the 
distribution of channel depth, providing an extra form of energy distri- 
bution not available in straighter reaches. This form of increased 
entropy is assumed to be significant in the stabilization of the 
meandering form, thus accounting for the meandering behavior of rivers 


(Leopold and Langbein, 1962). 
Dee PlantOnumo eres 


a) General 

Mean angular deviation and variance of angular deviations were 
derived for all reaches of all rivers studied. Table 4.03 presents 
first and second order moments for all the reaches. It is interesting 


to note that most reaches and most rivers, when regarded as separate 
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whole entities, show mean angular deviations approaching zero (marginally 
positive) when relative curvature, using both positive and negative values 
according to change in flow direction (derivation explained in detail 
in figure 3.14), is considered. At the sampling interval used, 

(Ax = one tenth mile) non-zero values indicate a preferred’ direction 
of deviation from flow path, and so may represent either a low frequency 
meandering tendency in a much larger series, or a definite flow trend. 

Indication of a flow trend would suggest the presence of a large 
scale process operating somehow within the landscape to cause flow 
elements to prefer one deviation direction over others. Since only 
local or relative curvature was derived, this implies that a non- 
random or deterministic element acting within or somehow underlying the 
fluvial processes is pulling flow in a preferred direction. A very 
large sample of deviation angles would be necessary in order to eliminate 
low frequency meander possibilities as the genesis of an apparent trend 
in the data. 

Detrending data before analysis effectively removes the 
influence of such a trend on statistics obtained in the analysis. It is 
interesting, however, to note and describe trends (apart from during 
spectral analysis) as these might help in the definition of underlying 


physical processes active in the establishment of the series obtained for 


analysis. 


b) Beaver River 
The Beaver River shows the most pronounced deviation preference 
among the rivers studied. The series for the reach concerned shows a 


mean deviation of 2.7 degrees to the right. Table 4.03 shows that both 
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the straight reach and the upstream meanders support this trend, 

while the downstream section of meanders exhibits a near-zero mean 
deviation. The total reach of river included in the study was rather 
small (only 18.5 miles of flow path), and it has been mentioned that the 
meanders recorded in the study are misfit in a large valley that meanders 
at a very low frequency relative to the oscillations in the Beaver's 
planform geometry. The study reach constitutes a small portion of a 
large concave-downward bend in the main valley, accounting for the 
apparent right trend in the series. 

Meanders on the Beaver River are marginally tighter in the 
upstream portion of the reach. Maps show that upstream from the study 
reach, the meanders tend to assume longer wavelengths relative to their 
amplitudes, while becoming more sharply curved in the downstream 
direction, shown by an increase in the mean angular deviation. 

In laboratory simulations of the intiation and development of 
river meanders, meanders begin to form at the furthest point downstream 
and the sinuous pattern moves slowly in the upstream direction (Friedkin, 
1945). Many scroll bars along the meanders, as shown in Figure 3.06, 
give evidence of the history of downstream movement of meanders. These 
have not been dated though, and it is not possible to interpret the rate 
of meander development from available data. Dating of scroll bars and 
data on much larger portions of the Beaver would be of assistance in 
the determination of patterns and processes operating on the Beaver's 


planform. 


c) The Kootenay River: 


Planform data for the Kootenay River suggest an overall tendency 
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toward equilibrium of angular changes in planimetric geometry. The 
slightly positive means derived could indicate a marginal preference 
for left turns (southwest flow deviation), following the geology of 
the region. 

Bends tend to become sharper in the downstream direction, and 
oscillation appears most intense in the lower reaches studied, as 


evidenced in increasing variance. 


d) The South Nahanni River: 

The South Nahanni River also shows a tendency to balance reversals 
of curvature, as its mean deviation is also near zero. Angular statistics 
in Table 4.03 suggest the possibility of a slight preference for an Wad phe 
turns, or north-west flow. This reflects the general orientation of 
the landscape in macro-scale. 

Bends tend to become tighter in the downstream direction in all 


channel types of the South Nahanni series. 


e) Theoretical Implications 

The near-zero deviations found for the South Nahanni and Kootenay 
reaches suggest a tendency for rivers to keep angular turns in directional 
equilibrium, implying the ability of these rivers to balance chance or 
determined initial deviations somewhere later in the series through the 
establishment of counter-deviations in the flow-path. Both the Kootenay 
and South Nahanni series are relatively large samples of their respective 
rivers. The Beaver series is, in contrast, a smaller sample of its river, 
and it exhibits a trend not found in the other series. A preferred 
deviation direction is expressed in the sampled planform, but superim- 


position of the series on a low frequency oscillation explains the tendency. 
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D. INFORMATION CONTAINED IN THE SPECTRAL AND CROSS-SPECTRAL STATISTICS 


Although the analysis of the depth configuration and planform 
of rivers, based on the first two moments of each of the respective 
series, allows some insight into the meandering process, the insight 
provided is rather limited and is little more than descriptive. Processes 
are identified, when at all, only in very general terms. Definition of 
relationship of variables is crude. Sophistication of analysis should 
provide more valuable information and more detailed insight. 

Spectral and cross-spectral analysis are described in greater 
detail in Chapter 2. Appendix 2 details calculations involved in the 
computer program used in the deviation of the statistics used. The 
standard "Biomedical Computer Program'! BMD:02T, expanded to take up to 
3000 cases, was used in all spectral and cross-spectral computation. 

All calculations were conducted using a 10 %Z mton ratio. 

The creation of a third depth series was necessary in order for 
the cross-spectral analyses to be performed. The corrected series were 
sampled at intervals of one-tenth mile along the river path. In order 
to minimize error of estimate by this method, at each tenth mile, three 
points (two in Series B) were averaged to obtain each of the depth values 
used to construct the new depth series. Each new series then contained 
exactly as many points as its corresponding meander plan series. This is 
necessary in cross-spectral analysis. A possible alternative to this 
method, which would have maximized resolution in the depth series during 
the sampling process, would be interpolation of the meander series, along 
some best fit curve between measured points. Resolution with Ax = 0.1 


miles, as in this study, is limited to processes manifesting wavelengths 
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of greater than 0.2 miles.only. Smaller ‘detail cannot be examined 
here. Interpolation of meander data would allow retention of detail in 
the depth series and would then provide more accurate analysis of the 
correspondence between meander plan and depth of flow. The sampling 
interval allowed adequate resolution for the purposes of this thesis, 
and so finer sampling and interpolation possibilities have not yet been 
explored. In any work requiring greater resolution, finer sampling or 
interpolation would have to be used. The sampling interval used here 
will obviously not pick out high frequency undulations such as dunes, 
but it is adequate in the description of large-scale processes within 
the river. 

The use of finer sampling or interpolation raises another 
problem, however. Most computer programs available for spectral and 
cross-spectral analyses will take only a specified number of cases. 

In order to use the techniques effectively, a suitable combination of 
degree of resolution and adequate length of record must be achieved. 
In most instances, this will likely require the expansion of the computer 


programs used. 


1. General 

The resolution chosen for this study limits the phenomena that 
can actually be examined to the large-scale processes characteristic 
of fluvial behavior. General patterns of bed configuration fluctuations 
can be analyzed and cross-examined with planform geometry, but specific 
small patterns are bypassed. Dunes cannot be examined or described 
using the new series data, but the identification of pool and riffle 


sequences will be within the scope of the analyses. Resolution of 
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wavelength is limited to those greater than 0.2 miles for Ax sel to 
Oe Lo mic 

Coherence between depth and curvature spectra are generally 
quite low. The use of absolute curvature rather than relative curvature 
seems to make little difference in the coherency. estimates. Resolution 


problems may thus be present in the coherency analysis. 
2. The Beaver River 


a) Depth Spectra 

Figures 4.19, 4.20, 4.21, and 4.22 present spectral results for 
the Beaver River. Depth spectra for the Beaver River exhibit stable 
patterns characteristic of all individual reaches. In all cases, 
predominant frequencies were in the intermediate range, with wavelengths 
from % to & mile contributing the majority of power to the spectra. 
A tendency for shorter wavelengths, or oscillations in depth of higher 
frequency, to become more important in the straighter section was 
observed. Higher frequencies also tended to contribute a greater 
proportion of variance in the downstream section of meanders, indicating 
the maintenance of the predominant mode of wavelength in depth series 
with an additional faeeeaeed importance of oscillations of shorter wave- 
length. No obvious trend of wavelength increasing downstream emerged 
in the depth spectra, but the sample was small, and the flow was confined 
within a much larger, older valley. Figures 4.19, 4.20, and 4.21 allow 
comparison of the three reaches in their depth spectra. Since the 
natural series are actually continuous, and records used in the calculation 


of spectral statistics are discrete, with a sampling interval of 0.1 miles, 
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A) PLAN AND DEPTH SPECTRA. | 
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FIGURE 4:19 
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A)PLAN AND DEPTH SPECTRA, PLOTTED IN LOG SCALE 
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A) PLAN AND DEPTH SPECTRA, PLOTTED In LOG SCALE 
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peaks and all frequencies listed are band limited. 

Spectra of individual reaches resemble those expected for 
second-order autoregressive processes, showing a peak at intermediate 
frequencies. In contrast, the depth spectra for the entire series 
(Figure 4.22) appears slightly more as a first-order autoregressive 
process, with power concentrated in lower frequencies. The longer 
reach demonstrates the contribution of much longer wavelengths (those 
approaching infinite wavelength are shown as the most powerful) reflecting 
the influence of the older valley. 

The straight section shows greater contribution of power by 
higher frequencies, suggesting oscillations of very short wavelength. 
Dunes were present in the river at the time of the survey, but resolution 
of the data does not allow detail to enter this larger pattern. 

All reaches showed a tendency to favor depth wavelengths of 
10 to 15 river widths, in conformity with the findings of Leopold, 
Wolman, and Miller (1964). These represent oscillations of about 37 
and 5 river widths, well within the range suggested by Yalin (1971) 
and approximating the distribution of regions of bed erosion, as Yalin 
(1972) suggested. 

Major wavelengths demonstrated a tendency to approximate a 
distance equal in magnitude to a value 7A for the reach, where A 
is the average cross-sectional plane parameter (assuming a rectangular 
channel). A is the product of average recorded depth and mean width 


of the reach considered. Table 4.06 illustrates the property found. 


b) Planform Spectra 


The spectra of angular deviations showed a similar pattern 
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for meandering sections. Variance and power patterns are more pronounced 
for relative than for absolute curvature, and absolute curvature tends 

to produce wavelengths twice those for relative curvature, as expected. 
Absolute curvature has a smoothing effect on the information contained 

in the planform series, and so shows a pattern of power more evenly 
distributed across frequencies. 

Again, spectra tended to resemble that for the second order 
autoregressive process, with predominant peaks at intermediate frequencies. 
The plan spectra for the entire reach reflected this pattern. Figures 
erro. 4.20, 4e21, "and 4.22 also give the plan spectra for” the series. 

Wavelengths of the meander data show a tendency to occur both 
at approximately twice the downstream distance of the depth wavelengths 
of their respective reaches, and at wavelengths equal to the wavelength 
of that of their corresponding depths. Figures 4.19, and 4.21 
illustrate this property for both the meandering sections of the Beaver 
River. 

Again, wavelengths shown as predominant by the spectra tend to 
approximate values of NA for their reaches (refer to Table 4.06). 

Plan wavelengths range between 6 and 20 river widths in distance, or, 


between about 27 and 67 widths, as suggested by Yalin (1971). 


c) Coherency Spectra 

Spectra of coherence-squared shows that the best relationship 
between depth and plan series occur at wavelengths equal to those found 
to be dominant for the depth series, as well as at wavelengths twice 
the value favored by the depth. The relationship is most striking in 


the straight section, as shown in Figure 4,20, although these coherences 
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arent, statietically significant. 
Most coherence-squared estimates are fairly low, however 
(less than .34, critical value for 90% tolerance computed according to 
Panofsky and Brier, 1963), and only the extremely high frequencies 
with relative curvature showed a significant wavelength. These are of the 
order of maximum possible resolution, or 0.2 miles. This contradicts 
the findings of Chang (1969) concerning absolute curvature, but supports 
direct observational evidence that plan and depth changes correlate 
spatially at high frequencies. Figures 4.19, 4.20, 4.21 and 4.22 show 
spectra and coherency using relative curvature for the same reaches using 
absolute curvature. Since, in general, absolute curvature does not 
appear to improve coherency estimates with respect to significance, and 
since direction of flow does seem to matter to the river in its activity 
(tendency to equalize reversals of curvature), relative curvature data is 
the basis of most coherency observations for this study. 
Coherence-squared wavelengths show the relationship to their 
mean rectangular cross-sectional areas. These dominant wavelengths are 
equivalent to 10 and 20 river widths, and 7 river widths, ranging between 


2 and 20 river widths. 


d) Low Frequency Waves 

Predominance of any low frequency activity seems de-emphasized 
by this analysis of specific reaches. Importance of power concentrations 
at long wavelengths is more apparent in the spectra derived for both plan 
and depth series of the entire reach studied. Waves ranging from 


infinite length to 1.7 miles in length are given more (power) in these 
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pectra. Power concentrations decrease toward intermediate frequencies 

and resume important contributions only at higher frequencies (1.471 cycles 
per mile and greater) or shorter wavelengths (0.7 miles). The higher 
frequency fluctuations of the Beaver River are evident in the smaller 
reaches, but the large scale activity is not obvious without the 

inclusion of much more data. Low frequency waves become more precisely 


defined with increase in sample size. (See figures 4.19-4.22.) 


3, Kootenay River 


Ficures 4:23, 4.24, and 4.26 show spectral ?results for the 
Kootenay River. A larger sample of points was obtained for the Kootenay 
than for the Beaver River. The Kootenay sample also included a larger 
variety of channel types and exhibited much less regular meandering than 
did the Beaver River. Both planform and depth series show a more 
even distribution of power among frequencies, emphasizing the presence 
of a number of wavelengths, none of which sufficiently summarizes the 
nature of the Kootenay series. Consideration of many wavelengths is 
necessary to gain an understanding of processes involved in the formation 
of the patterns derived. Table 4.07 summarizes the spectral findings 


for the Kootenay. 


a) Depth Spectra 

Depth spectra for the Kootenay show an increase in wavelength 
downstream along the river. This pattern is less obvious from maps 
but shows up clearly in the spectra of the bed configuration. Figures 
4.23, 4.24, 4.25, and 4.26 show this tendency. Goncentrations in 
power are clear at several frequencies. Lowest frequencies contain the 


majority of spectral power evidenced in the series, and power increases, 
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though less rapidly, towards very high frequencies. Three scales of 
activity between lowest and highest frequencies seem important. Wave- 
lengths of 0.9 to 0.7 miles, 0.5 miles, and about 0.3 miles stand out 
as intermediate spectral peaks in all reaches considered, among which 
the wavelength 0.3 miles seems most "powerful". These define the 
irregular meandering behavior of the Kootenay. 

Predominant wavelengths occur at between about 5 and 20 river 
widths, although some longer wavelengths are found with less power. 
Wavelengths for entrenched sections are apparently longer than those in 
less contined sections. Table 4.07 lists major spectral statistics for 
the Kootenay River. 

A preference for coherence wavelengths to occur as multiples of is 
evident in the table. Values between wand 8 Tlware favored in all but a 
preference towards one of odd or even number multiples is not obvious. 

All depth series appear to assume the form of a second-order 
autoregressive process. Peaks at intermediate frequencies interrupt 
the curve. All reaches exhibit this pattern, which is also consistent 


with observations on the Beaver's spectra. 


b) Planform Spectra 


Spectra of angular deviations also show increasing wavelength 
downstream for the Kootenay. Again, several intermediate scales of 
activity contain the power of the series. Predominant wavelengths 
occur as 1.0 and .5, .6 and .3, and .2 miles, with .6 miles concentrating 
more power than any other wavelength. This compares favorably with 
spectral power indicated by the Beaver, as planform wavelengths are 


twice depth wavelengths for the Kootenay as well. The pool-and-riffle 


1453 


jue brad wit eo snare 
dati scone albeit ITs nt aml: 
old sigh swadt Lad epi team sorsoR astia ‘e “ 
SeEn2TO0N ay to tolveied its , 

vith OS hue ® shade Gpawded) t2 qU950 aftgumfaveg amend: ot 
»tswoq eco! div tone? oc6 edjanelovaw xsarol emos ( 


= 
1 


it deed) catt vagnol \Teapteqqn 995 Gaeliosg Gndonstine 301 


oveW 
rr * 
tol aplseidess bavi aoygh Setem ooerl TO,.4 side? .eaoksoes Ii ; vor 


a] 1d enigitinw 28 29d0 01 wmddynsihyvey apenas aay : eee “a 
i 


i sud [le at Werove! urea §S bnawh avowied asudaV .efded ody mi 3¢ a 
** 

-euolvdo ton 2i eulad ition yeahh ase. a0 bbe to suo ebtswod | be —_ 
veb1d-jatusa s 36 mio? sty smuagpr 9) +que eeisea doqab ie 
Squivetni taforeuped? syntbamietial Ja eainet seadend phe . 

nsiniends onla ef dottw .wpetieg abide aidifaxe weronat [IA art 


> 


sBvioeqa a'weavesd of3 aa naokdaannol 
‘ a 


Aigaslyvew gyieeeion! wide Onis amet shin teliépae Yo 1d 3 a 
. to eadeoa ot nfonttindal levevee co ftgA \yerededt adg te wat 
eliganiavew tisnlaebetl .20lsoe oft Ya amie ada at 
wrk 4a rivwonoo valim a, ee feater: c, bre vee ove he oe me 
Adds pasa ainsi ait shaiaaahial Sails ib can an war 


aks 


wy 


o°T 1°S 
Oe £°9 
0°S oreo 
6°? 0°6 
aca rae 
1G. © SAE Il 
0°8 O>SZ 
66°T 97°9 
o* T 9°S 
Mm ff M 
ay: V 
yoo 
A 


1c°0 
TE0 
49° 


ee. 


a 


OS” 
0O0*T 


cc” 


1 


oa 


ex 


Cay 


M 


Vv 


(bp) Yet 


ag 


eo 


007 1 


MAALY AVNALOON HHL wos 


lo 
© 
re 
st 
¢ 
a) 
© 
N 
¢ 


T°? Leg 82°O 
O°? 5) cae 
L°G pea aoe 
vies BL 0) ia 


Le Sap iat 67° 


has T°ST ‘hb wesaisdyj» "¢ 
per BAS 1c 
”°Z cry Oe 
0°s O%Ge 00°T wearysuMog °Z 
LT S29 07° 
ie aS oC: 
aie £*97 0071 Sprezq -] 


Vv Rd 


syqdeq 


SOILLSILVLS TVYLOSdS GHZTTVeEaNao 


Ril Rad St 


= 
fou we 
aes 3 

ae 


ee: 


| 8 


it) 


Vee 
¥ 9: a 79" 
¥ 6° wes? ¥ 7° 
¥ ¢* ws volt 
ue an 
VS a Vv €° 
Ye Y Cia ase 
eel 
ae Yt 
ve ¥ Z° ee 


; Ware SV ri 
V/ Goo 10] SV 7 ee 


Fee <ur <a 
\O (a) N 


|< 
S 
4 


(<4 
N 


<q 
<t 


(< 
~ 


S jy [etq0eds 


Oy 


\<a 
st 


{<a 
~~ 


I< 
oo 


<i 
Ov 


i< 
+ 


ira 


S,\/ Teaqoeds 
yadeap 


daALa 


ccc 


IIc 


L1Te 


°ay 0 “4ApTA 
ueow 


AVNALOOM FHL 


Silepueow 
Ie [nsatit 


G29 weetyzsdn “¢ 


Slopurow 


le[nsetit 
0°6 weoijzsumop °Z 
Ce Spresq. "f 
‘ay ‘yadop 
ue ow 


MOI dIHSNOLLVISY Vi 80°? PTIPL 


ey 


a 


toe 


rv 
‘ 
o~ 


LGiLoeAe 


a) 


os 


ori 


*¢ ¥ 


‘ay 


) 


*e 
“af 
“so 


-_ 
= 4% 


STY.8 ebsc 


< 


“Vv 
i ae 
shec: 


cp’ te 


sequence postulated by Langbein and Leopold (1966) seems to hold for this 
semi-entrenched river, including reaches that vary from braided to canyoned, 
as well as for the regularly meandering sand-bed Beaver. Figures 4.24, 
4.25, and 4.26 demonstrate this characteristic. Table 4.08 presents 

the TA-wavelength relationship for the Kootenay River. The relationship 

is not as strong for this river, but it is still noticeable. Planform 
wavelengths occur as 5.4 to 16 river widths. Figures 4.23 to 4.26 give 


the plan spectra for the Kootenay. 


c) Coherency 

Table 4.07 shows a strong relationship between predominant 
wavelengths and /lw, as wavelengths appear to vary between about 277 
widths and 87} widths as multiples of 77. Coherence squared estimates 
are best at .3 to .4 and .6 mile wavelengths, consistent with the 
univariate spectral suggestions for the reaches. However, none of the 
estimates are significant at 90%. Marginal increase in wavelength is 
noted for downstream reaches even in the coherency estimates. 

Coherency also assumes the pattern of the second order auto- 
regressive process. Figure 4.23 to 4.26 demonstrates the similarity. 
Intermediate wavelengths contribute substantially to the power of both 


plan and depth spectra and this is reflected in the coherency estimates. 


d) Low Frequency Waves 


Presence of low frequency waves, having wavelengths in the order 
of 2 to 3 miles (and greater) is shown to concentrate some, though less, 
power in the Kootenay series than in the Beaver series, for both depth 
and plan. The presence of low frequency activity is not emphasized in 


analyses of the whole study reach. These waves represent oscillations of 
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the Order 10 mniver widths. Strong seaglogic controls «may account for much 


of the low frequency activity. 


4. South Nahanni River 
Profiles were obtained during both the upstream and downstream 
navigation of the South Nahanni River. Navigators of the boats believed 
a thalweg survey in either direction too dangerous for records, equipment, 
and personnel, and so the "anti-thalweg" was obtained both for the 
upstream and downstream traverses. The anti-thalweg is the channel of 
least velocity in the river. It now seems a thalweg survey would have 
been possible, and plans for future work include this possibility. 
Series B was surveyed as boats advanced downstream, and so it 
was a much faster survey, incurring loss of detail relative to series 
A. Navigation problems in the collection of series B also render it the 
least stable, and so observations are based on series A for the most 
part. Upstream surveys proved the least difficult throughout the study. 
Table 4.04 shows that depth data differ for the upstream and 
downstream surveys. Change in river stage during the time period over 
which ae surveys were conducted accounts for some change in depth, 
but it seems apparent that the series differences are also due to the 
direction of travel on the river and location of boat across the river. 
Two dimensional bed mapping would solve this problem, but it is costly 
and time consuming and was beyond the scope of this study. 
Series A was chosen as the most representative sample as flow 
stage, boat speed, and maintenance of desired route were all most 
stable during the upstream traverse. Series A also contains information 


for reaches further upstream than included in the Series B record. The 
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upstream profile also allows better data resolution, and so it forms 


the basis iof discussion for analysis of the South Nahanni River. 


a) Depth Spectra 

Depth spectral statistics for the South Nahanni River support 
the trends observed for the Beaver and Kootenay as wavelengths increase 
downstream along the path of flow. Comparison of the major meandering 
sections, First Canyon, the section including Second and Third Canyon, 
and Fourth Canyon, corroborate this theory even though they are all 
deeply incised canyons. This again seems to extend the scope of the 
observation. The processes involved in meandering seem to create 
similarities in bed configuration in spite of the variety of bed materials 
and bank constraints encountered. This implies that part of the meandering 
form can be attributed directly to streamflow processes, regardless of 
bed and bank materials. 

The pattern established as concentration of power for the depth 
series follows the pattern of power concentrations in planform of twice 
the depth wavelength is also present in the South Nahanni data. The 
pattern occurs in both the canyoned and braided sections and is shown 
in Figures 4.27, 4.28, 4.29 and 4.30. 

Power concentrations for the depth series occur in low frequencies 
and diminish as wavelengths decrease. This implies that longer wave- 
lengths of bed undulations characterize the bed form better than shorter 
wavelengths, suggesting that generating processes act in the South 
Nahanni on a large scale. Lower frequencies contain more power for the 
South Nahanni than for either the Kootenay or Beaver Rivers. 


Table 4.09 depicts the relationship between wavelengths and [Tw 
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A) PLAN AND DEPTH SPECTRA, PLOTTED IN LOG SC 
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FIGURE 4.27 
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A) PLAN AND DEPTH SPECTRA, PLOTTED In LOG SCALE 
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for the ides Wavelengths again tend to approximate multiples of Tw 
between about 17 and 20 river widths, consistent with the findings of 
VYelim C1972) % 

The depth spectrum for both the Second-Third Canyons and First 
Canyon best resembles a second-order autoregressive process, while the 
spectrum for the braided sections resemble more a first-order auto- 


regressive scheme. 


b)e’ Plantorm Spectra 

The tendency for plan wavelength to increase downstream is 
apparent in the South Nahanni River. Plan seems related to bed confi- 
guration where bed undulations occur mainly at twice the frequency of 
plan oscillations. Table 4.09 shows the relationship found. Power at 
similar frequencies also shows a relationship between plan and depth. 


In general, variance spectra of the South Nahanni's main reaches 


show power to be concentrated in both low and high frequencies. Occasional 


peaks in intermediate frequencies also concentrate power, and these are 
echoed in corresponding depth series spectra. Concentrations occur 
throughout the plan spectra at frequencies less than .385 cycles per 

mile and greater than 4.231 cycles per mile, representing wavelengths 
greater than 2.59 miles and less than .24 miles. Preferred intermediate 
frequencies include 2.308 cycles per mile (A= .43 miles) and 3.590 to 
3.718 cycles per mile (N= ,28 to .27 miles). The Limit of recolucion 
according to the sampling interval chosen is 0.1 mile, and so the analysis 
of wavelengths smaller than 0.2 miles was not possible from this data. 


Plan spectra show wavelengths of about 2 to 40 river widths to 
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occur in, the. spectra. 

Plan spectra typically appear as second-order autoregressive 
processes, with major power contributions within intermediate frequencies. 
Figures 4.27, 4.28, 4.29, and 4.30 show plan spectra for the Nahanni 


reaches. 


c) Coherency 

Low frequency coherence squares were generally high relative to 
those of other frequencies. Other preferred frequencies occurred as 
intermediate bands, and power in very high frequencies was small. 

Table 4.08 shows coherency wavelengths. 

A tendency for wavelengths to occur as multiples of w is 
evident. First Canyon, the most regularly meandering section, shows the 
relationship best. The coherency spectrum of First Canyon is presented 
in Figure 4.28. 

Activity appears concentrated in wavelengths of between 2 and 
20 river widths. Coherency is generally low using relative curvature 
data, and only a wavelength of .35 miles is significant at 90% tolerance 
limits. The use of absolute curvature improves the coherency somewhat, 
but still only three significant wavelengths, at Lacs cand 2 Cems 
are found. These represent distances of 3.4 and 1.1 river widths. 

Coherency spectra are multi-peaked but tend to resemble second- 


order autoregressive processes. 


5. General 
The spectra of all rivers studied generally tend to approximate 


second-order autoregressive processes, characterized by predominance of 
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intermediate frequencies within the data. This observation corroborates 
the findings of Ashida and Tanaka (1967), Nordin and Algert (1966), and 
Hino (1968) who proposed Markov second-order linear models for the 
Occurrence of sand waves and bed roughness. 

Wavelengths generally vary between about 2 and 20 river widths, 
consistent with Yalin's findings (1972). 

Pool-and-riffle sequences appear in the spectra for all types 
of reaches. The relationship, where plan wavelengths are twice depth 
wavelengths, is strongest in meandering and straight sections, but also 


occurs in braided channel. 
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CHAPTER V - CONCLUSIONS 


A. THEORETICAL CONTRIBUTIONS 

Classical methods of describing and explaining fluvial meandering 
were handicapped by subjectivity, partially due their lesser degree of 
technological and technical sophistication. Until recently, large 
amounts of data created more problems than they solved. The introduction 
of new techniques useful for both the collection and subsequent analysis 
of great quantities of data should prove beneficial to the science of 
fluvial geomorphology. 

The introduction of spectral and cross-spectral analytic techniques 
into fluvial geomorphology began in 1965 with Speight's work on the 
Angabunga River, but use of the technique as a descriptive statistic 
for fluvial systems seems to have advanced little in the past eight 
years. The technique allows examination of large amounts of data, 
but long series have not been available until recently. Although it 
is possible to record, continuously, depths along the path of flow of 
rivers, this possibility has not been widely exploited in fluvial 
geomorphology, prior to this study. 

A problem arises, however, in comparing spectral studies of 
meandering to the more traditional studies in the field. Meander 
wavelength was classically defined as down-valley distance. Spectral 
and cross-spectral analytic techniques require, for analysis, only 


single-valued functions. Rivers in nature can seldom be assigned a 
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baseline or origin allowing them this property in Cartesian space, 

and so, for spectral analytic studies, series are constructed down-path 
rather than down-valley. Resulting "wavelengths" represent down-path 
characteristics that cannot be compared with those found by traditional 
methods. More precisely, "path length" wavelength, as here defined, 
varies with wave amplitude, as defined by down-valley Cartesian co-ord- 
inates. The transformation between these two wavelengths is obviously 
the sinuosity ratio, but this can only be non-arbitrary (and objective) 
when meander wavelength is known. Thus, when multiple wavelengths 
co-exist, the sinuosity ratio is arbitrary (subjectively derived), and 
the transformation cannot exist (without serious bias causing its 
invalidation). 

Wavelengths derived in variance spectrum analysis, however, 
relate more directly the flow and its effects, as effects are measured 
along the path the river chooses in sculpturing its fluvial landscape. 
"Typical wavelengths, then, are no longer described, but entire 
stretches of flow path are analyzed as single continuous series. The 
occurrence of all possible wavelengths within the data is allowed, 
and even assumed, in spectral analysis, providing for the description 
of all possible scales of activity, above a size equal to twice the 
sampling interval, within the series. Relative contributions of 
different scales of activity are measured and defined. 

Intentions underlying the research conducted were focused on 
the collection and analysis of long, continuous series, sampled at 
discrete intervals, of plan and depth information for a variety of 
rivers. The actual study encompassed sand and gravel bed materials; 


large and small discharges; incised and freely meandering reaches; 
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braided, meandering, and straight channels; tranquil and turbulent flows; 
and glaciated and non-glaciated terrains. The research itself enconm- 
passed the employment, and to some extent, development, of both data 
collection and data analysis techniques. Some conclusions may be 


drawn: 


1. It is possible, and indeed desirable, to collect and analyze long 
records of plan and depth data, in order to attempt the discovery 


of "meaning" within the data. 


2. Boat speed varies along the river during the process of depth 
recording, and the resulting errors must be corrected. Systematic 
removal of points, in proportion to speed of recording, provides 
adequate and satisfactory correction for this problem. Detail 


can be retained if sampling intervals are sufficiently small. 


3. Digitization by machine (Digicon) of depth chart records, and 
subsequent storage of data on magnetic tape, is an extremely fast, 
accurate, and economic way of preparing data for analysis. In this 
study, about eighty-seven feet of depth chart, broken into 203 blocks 
for subsequent length correction, were diritized with, x — .02 duchess, 
so that over 52,000 (x,y) co-ordinates were stored, Cor a eotal 
cost, including Cal-comp plotting of the stored data, of about 


eighty dollars ($80.). 


4. Self-stationarity tests can be easily applied to data and provide more 
ready communication of information about the series for persons not 


involved in the research. Comparability of research is thus facilitated. 
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No previous researchers using spectral techniques for meander or 


depth forms seem to have tested for stationarity. 


Observations made with respect to the data collected using simple 
first and second moment statistics, and using variance spectrum analysis. 
Some theories (explained earlier in the review of literature) were 
tested. Reinforcement of some previous hypotheses via sophisticated 
statistical analysis, may be of some value to the discipline. Some 


conclusions may be drawn: 


5. Simple first and second moments may be useful in the detection of 
broad patterns of activity within the data. Underlying processes, 


such as geologic control, may be suggested by these statistics. 


6. Fluvial planform series tend to establish an equilibrium in angular 
deviation from flow path. Mean angular deviation tends strongly 
towards zero (0), or complete balance, on rivers with regular 
meandering, irregularly meandering, straight, incised, and even 


braided channels. 


7. Radius of curvature, based on circular paths, is a versatile 
statistic for the description of meander planform. It relates 
planform to flow acceleration in bends, and provides a non-dimensional 
plan parameter when divided by river width. The parameter R/w shows 
a strong relationship with the value l/f such that R/fiw appears to 
vary within narrow limits. Derivation of w makes difficult, in 


this study, precise definition of the relation, but the range from 
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12. 


1.1 to 2.3 seems to describe the fluctuations of R/yw in this 


study. 


Both plan and depth spectra indicate a tendency for down-path 
wavelength to increase in the downstream direction. Coherency 
estimates support this observation. The trend is not sufficient, 


however, to make the series non-self-stationary. 


Both plan and depth spectra indicate a tendency for most pronounced 
down- path wavelengths to lie within the range of about 6-20 river 


widths. 


Coherence wavelengths for the Beaver tend to occur as multiples of 
the product 7w (#times the river width). This is consistent with 


the findings of Church (1970) and the predictions of Yalin (1971). 


The Beaver River demonstrated a tendency for down-path wavelengths 
to gapproximate; 2 value 7A, where A is a mean cross-sectional area 
parameter based on the assumption of a rectangular channel. A is 
calculated as the product of the mean depth and mean width for each 
reach considered. This finding may be fortuitous, since the 


relationship was not as strong for the other rivers studied. 


Spectra of plan and depth and coherence tend to resemble that 
for a second-order autoregressive process, concentrating power in 
intermediate frequencies and decreasing in power towards high 
frequencies. These observations are consistent with the findings 


of Ashida and Tanaka (1967), Hino (1968), and Nordin and Algert (1968). 
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13. Down-path wavelengths for plan data tend to occur at values twice 
that for depth spectra, indicating the persistence of pool-and- 
riffle sequences in all manner of channels and bed material. The 
sequence is found for braided, straight, meandering (irregularly 
and regularly) and incised channels. Coherency spectra support 


this observation. 


14. Coherence is generally low between depth and relative plan data, 
although some significant coherences occur at greater than 90% 
probability. In general, coherence is not greater when absolute 


curvature data are used rather than relative curvature data. 


15. Similarities in bed and plan configuration occur for rivers of 
various bed and bank materials, indicating a more direct relation- 


ship of form to flow than to bed and bank media. 


16. Spectral and cross-spectral analytic techniques increase objectivity 
while maximizing quantities of data analyzed, and results can be 
related both to other variables and to the map. Wavelength is 
defined as down-path distance, measured along the path of flow. 
Coherence SEE, can provide insight into the relations between 


form and process. 


B. RECOMMENDATIONS FOR FUTURE RESEARCH 

The research conducted in the completion of this thesis was 
in many ways new. New methods of data collection and preparation were 
attempted, and relatively new methods were employed in the analysis of 


the data. Many problems were encountered in the endeavour, and though 
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some were solved in the process of this work, others remain unsolved. 
Ideas encountered are offered below as recommendations towards research 


into the broader problem--the mystery of fluvial meandering. 


1. Improvement of Type of Data Collected 


Accurate measurement of river widths recorded in the field at 
the time of the depth survey would be useful in the examination of 
precise relations between hydraulic variables. Widths were obtained 
from maps and air photographs for this study, and the limits on their 
accuracy restricted the use of some ideas that seemed promising. Widths 
could be measured for both bankfull and flood stages. 

Velocity measurements could be used in the anlysis of relation- 
ships between other variables of interest, and they could also be useful 
in data correction, accounting for boat speed relative to banks. 
Accurate measurement of velocity at a variety of depths and cross~-river 
positions would be valuable to an understanding of many physical 
processes acting within the river to form the depths and meanders 
observed. 

Sediment samples would be of further use in the definition of 
processes operating in the river, and they would facilitate comparison 
of reasearch. A knowledge of sediment grain size would broaden the 
scope of the work. 

Depth and plan data should be recorded so that they correspond 
as precisely as possible. Where possible, thalweg data should be 
obtained first. 

Slope data would also be useful in the examination of fluvial 


meandering processes. Information on both water surface slope and bed 
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slope would be extremely valuable. 

Data gathered at bank altering stages of flow would be valuable 
in solution of some of the problems involved in meandering. Ready 
access to boats and recording equipment is necessary, however, for this 


to be accomplished. 


2.yeResolutionvof Data 

Where digitization of data is little problem, as much data as 
possible should be gathered. Recording charts should be run quickly 
to maximize information obtained. 

In order to preserve maximum resolution of depths in cross- 
spectral analysis, planform data must be expanded. Interpolation is 
suggested when analysis is to focus on small-scale processes within the 
river, although data removal proved most satisfactory toward the exami- 
nation of large-scale activity in this study. 

Optimally, two-dimensional mapping of bedforms would eliminate 
many problems of the depth survey while expanding scope for data 
analysis tremendously. Instrumentation exists for continuous mapping, 
but is not readily available, owing to its expense. In lieu of this 
degree of sophistication, a number of channel cross-sections could be 
obtained to improve utility of data and provide some of the information 
sought. This step of data collection requires time and labor, and thus 


incurs additional expense. 


32) | Accuracyvof Absolute: Location 


Absolute location must be maintained as accurately as possible 


during the depth survey. Air photos and. maps are not always sufficient 


for this purpose. Regular recording of compass bearing and time elapsed 
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since beginning of survey, with respect to chart position, could help 


to overcome this problem. 


4. Data Analysis 


Programs for spectral and cross-~spectral analysis must be able 
to handle vast amounts of data if interpolation methods for retention 
of resolution succeed. Dr. M. Church at the University of British 
Columbia (M. Church, pers. comm.) is presently completing work on a 
program designed to handle up to 100,000 data points in the eventual 


solution of the problem of meandering. 


5. Expansion of Breadth of Research 


Data must eventually be gathered and analyzed for rivers of 
all sizes and channel forms. Co-operation among workers in the field 
is essential for this to happen. An exciting possibility, based on the 
work of Yalin (1971), is the concept of obtaining two-dimensional 
spectra of bedforms, and relating these, via 2-D cross~spectra, to two 
dimensional water velocity variations. If the hypothesis that differing 
scales of turbulence produce observed bedforms is correct, then phase 
spectra, transfer functions, and coherence spectra should provide 
elucidation of these processes, and the writer plans to do further 
research in this area. 

Applicability of results obtained will be known more precisely 


as these suggestions are carried out, and prediction of form and process 


will become more possible. 
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APPENDIX 1 


STATISTICS FOR DEPTH SERIES CORRECTION 


(5) COMPUTATIONS IN CORRECTION 


~ Standardization of Series 


Step l. Calculate mean block speed, V, for each block, as: 


y = distance travelled 
time elapsed 


= Map oretance | © ig UI aes 
chart length L/L T 
chart speed r 


Step 2. Form ratios of each mean block speed, v;, to fastest block 


1c 


speed found, Veen, as: 
vi 
ee ces 
Vmax 
Steu >. Remove, via systematic sampling, 100(l-r) per cent of the 


points in each block. 
Step 4. Join all points remaining, in sequence, to form the 


"standardized" depth series. 
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APPENDIX 2 
COMPUTATIONAL PROCEDURE USED IN SPECTRAL 
AND CROSS~SPECTRAL ANALYSES OF DATA 


(source: BMD:02T program) 


COMPUTATIONAL PROCEDURE 


1. Symbols Used 
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iralal ; é ; ; 

i value of discrete, equi-spaced time series X(t) 
Agee" : : } 

it value of time series X(t) after subtracting mean 
th , ; : 

i value of time series Y(t) after subtracting mean 
i value of time series after prewhitening 

Lag 

Constant value (JCG}] < 1.0) provided by user 

Number of discrete data points 

Maximum lag 


Constant time interval 


: . Pp 
Autocovariance of series X at lag p, (2° 


At 
ye ng? , 
Autocovariance of series X after detrending at lag p 


Power spectral estimate of series X at frequency 


jal 


se eee) 


Smoothed power spectral estimate of series X at 


frequency h7 
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Power spectral estimate of series X at frequency 


h@ after recoloring 
mAt 


Cross-covariance between series X and Y at lag p 


Cross-covariance between series X and Y after 
detrending at lag p 


Cross-power spectral estimate at frequency h# 
mAt 


h 7 
Cospectrum, the real part of pat frequency hi 
; At 


=) 


(h) 


Quadrature spectrum, the imaginary part of Poy 


at frequency h7 
mAt 


Smoothed cospectrum at frequency hI, 
mAt 


Smoothed quadrature spectrum at frequency he 
mAt 


Amplitude of cross-spectrum at frequency h7_ 
mAt 


Phase of cross-spectrum at frequency h@/ 
mAt 


Transfer function at frequency h/7 
mAt 


Amplitude of transfer at frequency h// 


mAt 


Coherence square at frequency h7 
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Procedure 
Step 1. The mean of the data is first calculated and subtracted 


from each value of —. 
n 
Rese! ; 
mec > Xi 
i=l 


Ae eye me i Ne ae eee, FT 


This is not done if the series is to be detrended. 


Step 2. (optional) One method of prewhitening the input data is 


by the following moving linear combination: 


C is a value supplied by the user ( |[C[{<1.0) 


Step 3. The autocovariance is then computed. 


G a 
y 1 
a8 ae = ee 
R, n-p Oe ee 0, 3 2 5 TH 
qui 
(p) 
Reo R, Cpt ) 


Step 4. (optional) If the series is to be detrended, a least 


squares fitting method is used 


cP? bi ne 


~§ + xivwhered = opiliew.., n=), where 
x x 


nh 
rig Weas x, (2i-nt1)/(C(n-1L)n(nt1))/6) and § = x - <(n-1)/2 
i=0 


where X is the mean. 
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Step 5. 


Step 6. 


DEED fs 


Step 8. 


ee res ) A 
P. = Sey cos nt, h = a Ie eee, TM 
1 0 <p<m 
where€ = 
Li 2ee se Oy a 
h 
ites ee ee (a) = P (az. | 
mAt 


The raw estimates of the power spectrum are then smoothed 


by "hamming". 


(Oye (0) (1) 
spih = Eas 1) + .54p(h) oF ee 0 Zhen 
x x 
(m) 1 
sp™ = SAP m rm Aor ) 
x x 


A check sum is computed to check the computations of the 
estimates 


m= 1 
cuxsum = — Ee cept? + sp de oH] 
mAt 2 nel x 


(0 


This should equal R , the autocovariance at zero lag. 
x 


To compensate for the prewhitening in Step 2, the smoothed 


spectrum is recolored by 


i) 80D eF ae 
RSP. = oF A NG =e COS aie “At) 
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Similarly, Steps 1 to 8 can be applied to series Y. 
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Step 9. The cross-covariance is computed by the formula 


pe? eter Ps Sete: B= Oy Ly ey ees 


q= 
eS 
Crip)” oak . 
Ry ee Xotp ae Dia hy 2 eae gt 
q=l 


Step 10. (optional) In detrending the series 


COD ica Ge, aes 
My = Ry -p- OOP eres antl) La ers ot koe 


“and p are defined in Step 4. 


Step ll. The .crocse—spectrum, ee is given by 


Cy a 
iy tie Syyerti Oey 


The cospectrum is obtained from 


m 

Ce RG Roe at €p) (-p) ene 
= =: “i 2 

Coy oe Rey b oe COS = 


The quadrature spectrum is obtained from 
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xy Y pe0 P\ xy xy 
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Step 12. Both the cospectrum and quadrature spectrum are 


smoothed by “hamming'' as shown in Step 6. 


iat ae 
amir? - Boy F +f te the amplitude 
phase of the 
cross-spectrum 
are computed 


a (h) se, 
PHAS = Arg(SC,. + Say ) 
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The phase angle computed is the phase shift, of sernes. 


with respect to series X. 


Step 14. The transfer function from series XvboeYiisigiven by 
(h) 


Vee ea 

sie RV ok iPHAS 

Tyy pth) LNs e xy 
x 


The amplitude is given by 


h (h) (h 
xy xy x 
The phase shift from X to Y is the same as that of the 


cross-spectrum. Both the transfer function amplitudes 


from X to Y and Y to X are printed. 


Step 15. Finally the coherence square, between X and Y, is given as 


2 
cosy” = ( at) / se . gph? 
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